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Abstract. Physics of the Poynting-Robertson (P-R) effect is discussed and compared with the 
statements published in the past thirty years. Relativistically covariant formulation reveals the 
essence of the P-R effect and points out to nonphysical explanations in scientific papers and 
monographs. Although the final equation of motion 
m dv/dt = ( SA'Q' pr I c ) {(1 - v ■ e/c) e - v/c} 

has been usually correctly presented and used, its derivation and explanation of its essence is 
frequently incorrect. 

The relativistically covariant form of the equation of motion yields the P-R effect as an action 
of the radiation pressure force on a moving spherical body. No "P-R drag", as a particular 
relativistically covariant equation of motion, exists. Omission of the nonphysical term "P-R 
drag" excludes any confusion in the published definitions. 

The difference between the effects of solar electromagnetic and corpuscular (solar wind) 
radiation is stressed. The force acting on the particle due to the solar wind (the simple case of 
radial solar wind velocity is considered) is 
F sw = F 3W [ ( 1 - v ■ e/v sw ) e - x' v/v aw ] , 

where F sw is the force on the stationary particle, v sw is the heliocentric solar-wind speed, and, 
the value of x' depends on material properties of the particle (1 < x' < 3). 

Secular evolution of orbital elements is presented. Initial conditions are included. 

Key words, cosmic dust - radiation pressure - equation of motion - relativity theory - Poynting- 
Robertson effect - solar/stellar wind effect - celestial mechanics 
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1. Introduction 
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Motto: 

What is popular is not always right; 
what is right is not always popular. 



The action of electromagnetic radiation on a moving spherical body is called the Poynting- 
Robertson (P-R) effect. It is used in astrophysical modelling of orbital evolution of dust grains for 
many decades (e. g., Poynting 1903, Robertson 1937, Wyatt and Whipple 1950, Robertson and Noonan 
1968, Dohnanyi 1978, Burns et al. 1979, Kapisinsky 1984, Jackson and Zook 1989, Leinert and Grim 
1990, Gustafson 1994, Dcrmott et al. 1994, Reach et al. 1995, Murray and Dcrmott 1999, Woolfson 
2000, Danby 2003, Quinn 2005, Harwit 2006, Grim 2007, Sykes 2007, Kriigcl 2008, Wyatt 2009). The 
P-R effect holds when material of the particle is distributed in a spherically symmetric way. It is 
assumed that the spherical particle can be used as an approximation to real, arbitrarily shaped, par- 
ticle. Relativistically covariant equation of motion for arbitrarily shaped dust particle can be found in 
Klacka (2004, 2008b), where also thermal emission of the particle is taken into account (Mishchenko 
2001: the radiation pressure on arbitrarily shaped particles arises also from an anisotropy of thermal 
emission). Krauss and Wurm (2004) have published the first experimental evidence that nonspherical 
dust grains behave in a different way than the spherical ones. 

We will deal with the P-R effect in this paper. It has elapsed thirty years since the time of 
Dohnanyi's statement "Much confusion in the literature has existed on this subject and to some 
extent still does, unfortunately." (Dohnanyi 1978, p. 563). The situation has not changed and many 
"explanations" of the P-R effect have been offered during the past thirty years. It is important to have 
correct explanations of the effect, in disposal. 

The paper presents many statements on the Poynting-Robertson effect, which exist in the liter- 
ature. We have tried to show physically incorrect points of the statements and to put them into a 
correct way. The enormous amount of incorrect presentations probably confirms the ideas of Weisskopf 
(1973): "The pressure of fast publication is so great that people rush into print with hurriedly written 

papers and books that show little concern for careful formulation of ideas Essential ideas of physics 

are often lost ..." 

We have already tried to explain the physics of the P-R effect (see, e. g., Klacka 1992a, 1993b, 
2000b, 2001b). However, our comments and explanations have not met, in general, with a positive 
reaction and the incorrect statements and derivations are repeatedly cited and reproduced in many 
other scientific papers and textbooks. This avalanche propagation is the reason why we have decided 
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to write a paper presenting comments on various derivations and explanations. This will lead to a 
deeper and more complete understanding of the effect. 

Our paper is based on the results of Klacka (2008a, 2008b). At first, we summarize the most 
important results on the P-R effect. Then, we discuss many "explanations" of the P-R effect, which 
have been presented in the literature during the past thirty years. The scientists should concentrate on 
Sec. 6.1, where comments on the most frequently cited paper by Burns et al. (1979) are given: not only 
the important parts of Klacka (1992a - Sec. 2.5) are repeated, but the relativistic results of Klacka 
(2008a, 2008b) are used. It is clear that the statement "the interpretation of the low-order velocity 
terms in classical (nonrelativistic) physics are still debated" (Gustafson 1994, p. 565) is already not 
relevant, now. Really, even the simplest case of Poynting (1903) and Robertson (1937) is correctly 
understood on the basis of relativistic approach. Moreover, the dimcnsionlcss efficiency factor Q' vrl 
Q' ext are obtained from Mie's solution of Maxwell's equations and it cannot be said that Maxwell's 
equations belong to classical (nonrelativistic) physics. These relativistic results are used also in other 
parts of the paper; e.g., the case of the aberration of light is discussed in detail in Sec. 6.2. 

To be sure that the physical results are correct, one has to consider relativistically covariant equa- 
tion of motion. Klacka (2008a, 2008b) explicitly showed the advantage of relativistically covariant 
formulations: the application of this approach enabled to found that even the statement of Poynting 
(1903) and Robertson (1937) (and other explanations up to the year 2008) is incorrect: as for the out- 
going radiation, the particle's "own radiation outwards being equal in all directions has zero resultant 
pressure" (Poynting 1903) in the particle's own frame of reference, or, "the process of absorption and 
re-emission produces no net force on a particle when one chooses to work with a stationary frame 
referred to the particle" (Wyatt and Whipple 1950). 

The importance of relativistically covariant equations of motion is stressed, e. g., by Feynman and 
Mould. Feynman states: "The laws of physics must be such that after a Lorentz transformation, the 
new form of the laws looks just like the old form." (Feynman et al. 2006, volume II, p. 25-1), or, 
"There is, however, another reason for writing our equations this way. It has been discovered - after 
Einstein guessed that it might be so - that all of the laws of physics arc invariant under the Lorentz 
transformation. That is the principle of relativity. Therefore, if we invent a notation which shows 
immediately when a law is written down whether it is invariant or not, we can be sure that in trying 
to make new theories we will write only equations which are consistent with the principle of relativity. 
... the interesting physical thing is that every law of physics must have this same invariance under the 
same transformation. ... It is because the principle of relativity is a fact of nature that in the notation 
of four-dimensional vectors the equations of the world will look simple." (Feynman et al. 2006, volume 
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II, p. 25-11). Feynman has motivated a reader to the importance of using tensors (or, as a special case, 
four-vectors) in formulation of equation of motion. In order to be the reader more familiar with the 
relevance of this access, we reproduce the corresponding explanation from Mould (2002, p. 123): "It 
is a goal of relativity to formulate the fundamental laws of nature in covariant form; that is, in a form 
that is independent of coordinates. A four-vector equation, for instance, is covariant in this sense. The 
particular values of a four-vector may very well depend on coordinates, but a four-vector equation will 
have the same form in all coordinate systems. Tensor equations have this same property. They have 
the same form in all coordinate systems even though their specific component values differ from one 
coordinate system to another. Scalar functions have a form invariance as well as value invariance." 

Sections 2-5 present the most relevant results on the physics of the P-R effect. Sec. 6 discusses 
various "explanations" of the P-R effect, published during the past 30 years. Sec. 7 offers two simple 
correct explanations of the P-R effect. Finally, Sec. 8 summarizes the relevant equations for the 
secular evolution of orbital elements of a spherical body under the action of central star's gravity 
and electromagnetic radiation. Discussions presented in the paper can help in understanding the P-R 
effect. The paper should ensure that only correct explanations will be published in future. 



2. Covariant formulations 

Relativistically covariant formulations are useful in understanding of the interaction between a particle 
and an incident electromagnetic radiation. Spherically symmetric mass distribution within the particle 
is tacitly assumed. 

Interaction of the spherical particle with the incident electromagnetic radiation is known as the 
Poynting-Robertson effect. The emission component of the radiation force, due to thermal emission, 
is zero, i. e. F' e = 0. The conditions for the incoming and outgoing radiation are (Klacka 2008b): 



dp? _ w*SA'Q> ext 



dr c 
dvLt = w 2 SA'Q ' ext ( »" 
dr c 



[q'^ + (1-QO &"} , 



eg - 5^ , (i) 

*°t ext 

where the left-hand sides are four-momenta, e.g. p^ ut = (E out /c ; p out ), c is the speed of light in 
vacuum, S is the flux density of radiation energy (energy flow through unit area perpendicular to 
the ray per unit time) , A' is geometrical cross-section of the spherical particle, Q' ext is dimensionless 
extinction efficiency factor given by optical properties of the particle and averaged over the stellar 
spectrum, Q' is dimensionless efficiency factor for radiation pressure integrated over stellar spectrum 
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and calculated for radial direction (as for the dimcnsionless factors of effect ivity of radiation pressure 
and extinction, see Mie 1908, or, e. g., Klacka and Kocifaj 2007, or, Sec. 4.5 in Bohren and Huffman 
1983) and 4-vectors u M (four- velocity), 6 M and the quantity w are given by the following equations: 

w M = (7 c ; 7 v) , 
1 e 

w ' w 



w = 7 



where v is particle's velocity in the reference frame of the star and e is unit position vector of the 
particle with respect to the star; 7 is the Lorentz factor. As was pointed out by Klacka (2008b), Eq. 
(1) differs from the geometrical approach presented in the literature since the time of Poynting and 
Robertson. 

Equations for dp^ ut /dr and dp^ n /dr yield the following covariant relation between the outgoing 
and the incoming radiation: 

^ = (1 - Q>) *L + w2SA 'Q'ext q, * 

Q' = ■ (3) 

ext 

Eq. (3) can be considered as a fundamental equation for the validity of the P-R effect. 
2.1. Incoming radiation and particle's acceleration 

The action of the incoming radiation influences particle's motion according to the following equation 
of motion: 

dp» dpi 



in 



dr dr 

where the right-hand side is given by Eq. (1) and 



(4) 



p^ = mu^ (5) 

is a four-momentum of the particle. Inserting Eqs. (1) and (5) into Eq. (4), we obtain 

du^ dm u w 2 SA'Q' ext 

m — 1 — - — = -?L_ext_ y _ (e) 

dr dr c 

Multiplication of Eq. (6) by u u yields, on the basis of u M du^/dr = 0, u n = c 2 , u u = c: 
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Inserting Eq. (7) into Eq. (6), we obtain four-acceleration 

du»\ _w*SA'Q> ext ( b ,_^\ (8) 



dr J m mc \ c . 

The result represented by Eq. (7) can be understood in the sense that dE in /dr = w 2 SA'Q' ext 
corresponds to the energy of the incoming radiation interacting with the particle [dEi n /dr = c? 
(dm/ 'dr)i n ]. The real system consists of radiation and the particle and the system cannot be divided 
into the individual components. Thus, the result presented by Eq. (7) does not correspond to the real 
change of the particle mass, but only to a virtual change. No direct real experiment can verify the 
validity of Eq. (7). Indirect evidence of Eq. (7) is the validity of the result < Q' pr /Q' ext < 2 and 
experimental confirmation of Maxwell's equations and relativity theory. 

It may seem surprising that (dm/dr) in depends on the particle velocity and A'Q' ext , only. E.g., 
perfectly reflecting and absorbing spherical surfaces, within geometrical optics approximation, yield 
the same result for (dm/dT) in . 

2.2. Outgoing radiation and particle's acceleration 

The action of the outgoing radiation influences the motion of the particle according to the following 
equation of motion: 

dr ~ dr ' 1 ' 

where the right-hand side is given by Eq. (1). Inserting Eqs. (1) and (5) into Eq. (9), we obtain 



in 



(XT (XT C I C J 



Q'^p-- (io) 

V ext 

Multiplication of Eq. (10) by yields, on the basis of u a u M = c 2 , u a du^/dr = 0, and, u a 6 M = c: 

dm\ w*SA'Q' ext ^ (n) 



dT / out 

Inserting of Eq. (11) into Eq. (10), we obtain four-acceleration 

(duP\ _ w^SA'Q' ext (l _ Ql) ( b ,_u^ 
\ dT ) out mc \ c 

Q' = -fir 1 - ■ ("J 

V ext 



2.3. Equation of motion 
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Equation of motion of the particle is given by the condition 

(13) 



dp» drf n dp^ u 



dr dr dr 
Eqs. (1) and (13) yield 



Q> 6" - - 



dp^ = w 2 SA'Q> ext 
dr c ~* \~ c 

Q' = ■ ( 14 ) 

W ext 

Insertion of Eq. (5) into the left-hand side of Eq. (14) yields dp^/dr — m{du ti /dr) + (dm/dr)^. 
Using this relation and subsequent multiplication of Eq. (14) by u^, one obtains dm/dr = 0, if also 

(du^/dr) — and — c 2 , u^b^ — c are used. The relation dm/dr = states that the mass 

of the particle is conserved. Thus, Eq. (14) may be rewritten to the form 

du^_ = w*S A' Q' ext f 
dr m c \ 

Q' = jsr- ■ (is) 

W ext 

3. The fundamental condition 

In order to be able to correctly understand physics of the P-R effect, we have to realize that the 
fundamental condition for the validity of the P-R effect exists (Eq. 3). The essence of the P-R effect 
is that nonradial components of the particle's radiation pressure and components of the thermal force 
are zero (in particle's frame of reference). Thus, the P-R effect holds for spherical particles. 

The fundamental condition states: the total momentum (per unit time) of the outgoing radiation 
is nonzero, in general, and 

F' - F' 

^outgoing ^incoming 5 



P outgoing — P incoming 

Q' 



Q pr 



ext 

F' e = 0, (16) 

in the proper frame of reference of the particle. The condition is obtained from its covariant form 
presented by Eq. (3). The equation states that the P-R effect holds if the total momentum (per unit 
time) of the outgoing radiation, integrated over the whole space angle, is in the direction of the incident 
radiation, in the particle's frame of reference. 
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It is generally believed that the special case Q' = 1 corresponds to the perfectly absorbing sphere, 
in geometrical optics approximation. The equation of motion for this special case was discussed by 
Poynting (1903) and finally written by Robertson (1937). It is believed that this special case yields 
i) p' 'outgoing = dp' out / dr = (sec Eqs. 1 or 3) and this corresponds to the hypothesis of "isotropic 
reemission of the radiation", and, ii), the case Q' = 1 is also produced by specular reflection on the 
spherical particle (within the geometrical optics approximation). However, the physical condition for 
the two cases, perfect absorption and specular reflection, is Q' = 1/2, in reality (Klacka 2008a, 2008b). 



4. Approximation of the equation of motion 

It is sufficient to approximate space component of the covariant equation of motion (and also acceler- 
ation of the particle) to the first order in v/c. Eq. (14) yields the following relation between its time 
and space components: 

dp = v dp 
dr c dr 

For this reason, we have to do an approximation of the time component to the second order in v/c. 

At first we will present approximations for the incident and outgoing radiation. Eqs. (1) yield the 

following approximation (the results are exact, accurate to all orders in v/c): 

dE ln „ ., -. ( _ v ■ e> 
~dt 
dp lr 



— S A Q' ext ^ 
' ^ 1 



1 



SA'Q', 



c 

v ■ e 

c 



2 / 



dt c 

The results for the outgoing radiation (to the second order in v/c in energy) are: 

v ■ e /v ■ e N 



(18) 



dE out 
dt 



— S A' Q' ext 



(- 



dpout = SA'Q' ext 
dt 



Q' 



C 
pr 



(19) 



Q'ext 

where t is time measured in the stationary reference frame, i.e. the reference frame in which the 
particle moves with the velocity v. Acceleration of the particle, due to the incoming and outgoing 
radiation, we get from Eqs. (8) and (12) 

v ■ 



fdv 




SA'Q' ext 


[Tt 


^ in 


mc 


'dv\ 




SA'Q' 


.dt) 


out 


mc 




Q 1 = 


Q 'pr 
Q'ext 



ex I 



V 

C 

v ■ e 



(20) 
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Covariant equation of motion of the particle yields 

dp = SA'Q' ext r/ _ u_ex _ tn 
di c IV c / c J ' 

Q' = |^ (21) 

V ext 

and the acceleration of the particle is 



dv SA'Q' 



ex t 



dt mc 

Q' = %^- (22) 



ext 

5. Radiation pressure and the Poynting- Robertson drag 

The P-R effect is characterized by the following equation of motion to the first order in v/c: 
dv GM r / v ■ e\ 

dt U 1 - — ) e -J ' 

4nGMmc ' 

where G is a gravitational constant, L and M are luminosity and mass of the central star around 
which the particle moves, r is a distance of the particle from the star and R is a radius of the spherical 
particle. 

Gustafson (1994, p. 566) writes: "The velocity-independent radial term representing force [acceler- 
ation] due to radiation pressure... - although not a pressure, it is often referred to simply as 'radiation 
pressure' " : 

dv\ n GM 

0-Z-e. (24) 



^ / pressure 

According to a definition (see, e.g. Gustafson 1994, p. 566, or, Dohnanyi 1978, pp. 562-565), the 
velocity-dependent part of Eq. (23) is the P-R drag: 

uL B „ = -"-{(-) « + c)- < 25 » 

We can decompose the velocity vector v of the particle into the radial and transversal parts 
v = vr e + vt &t = r e + r8 ex , (26) 

where the unit vector ej is normal to e in the orbital plane. Substituting Eq. (26) into Eq. (23) yields 

dv n GM f / „ f\ r() \ 

= P^T il 1 ~ 2 7 e - -Ct} ■ (27) 



c I c 



10 



Klacka et al.: Poynting- Robertson effect 



On the basis of Eqs. (25) and (26) we can write for the P-R drag 



(—\ - (3 — ( 2 - e + — e \ 

\ dt J P-R drag \ C & C / 



(28) 



We have to stress that any decomposition to the "radiation pressure" (Eq. 24) and the "P-R drag" 
(Eqs. 25 and 28) is nonphysical. As follows from covariant formulation, the physical equation is Eq. 



Moreover, avoiding the usage of the "P-R drag" not only represents correct physics, but any 
confusion about the "P-R drag" is removed. This is also important since, besides definition represented 
by Eq. (25), also another "P-R drag" is used. It contains only the term proportional to — v/c in the 
equation of motion - the first term in brackets on the right-hand side of Eq. (25) is omitted (see, e.g., 
Wyatt 2009). One must bear in mind that the two different definitions of the P-R drag lead to two 
different results for secular changes of orbital elements [dap/dt = — (3 (GM/c) Ap-p drag / (ft/3 (1 — 
ep) 3 / 2 ), dep/dt = -13 (GM/c) E P - R dra g / (aj (1 -ep) 1 / 2 ), where A P _ R drag = 2 + 3e^ E P - R drag 
= (5/2) ep for the first type of the P-R drag and Ap-p drag = 2 + 2 e 2 3 , Ep-p drag = 2 ep for the 
second type of the P-R drag; ap and ep denote semimajor axis and eccentricity of the particle's orbit 
around the central star, see Eqs. 50-67 in Klacka 2004]. Unfortunately, another definition of the P-R 



see, e.g. Minato ct al. (2004). Really, confusion in the literature exists. And physics is ignored. 

The P-R effect is the effect of radiation pressure on a moving spherical (nonrotating) body. Thus, 
any separation between the radiation pressure and the P-R effect is also nonphysical (see, e.g., Leinert 
and Grim 1990; Mann 2009; see also Minato et al. 2004 - Sec. 3.2). 

6. P-R effect and statements presented in the literature 

As we have already mentioned at the beginning of the paper, much confusion exists in explanations of 
the P-R effect. As for the papers by Poynting (1903), Robertson (1937), Wyatt and Whipple (1950) 
and Burns et al. (1979), the physics was partially commented in Klacka (2008a, 2008b) and also in 
Klacka 1992a (Sec. 2.5). We will deal mainly with newer explanations presented in the literature, in 
this section. 

We have to bear in mind that the results obtained by Klacka (2008a, 2008b) are consistent with 
Maxwell's equations and relativity theory. Since relativity theory was motivated by better under- 
standing of Maxwell's equations, we consider Maxwell's equations (and their consequences) as an 
indispensable part of the relativity theory. Classical physics considers that interactions spread with 



(23). 



drag exists: 

(dv/dt) P -R dr. 



■ag — 



/3 (GM/r 2 ) {r6/c) e T , 
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an infinite speed. Thus, all the results presented in Klacka (2008a, 2008b) have to be considered as 
a consequence of the relativity theory. Any statement of the type "radiation forces - including the 
Poynting-Robertson drag - are fundamentally classical forces and are not produced by relativistic 
effects as commonly believed" (Burns et al. 1979, p. 8) is not physically correct. As an example, 
we can mention the calculation of Q' pr . It is based on Mie's solution of Maxwell's equations and it 
cannot be based on classical physics inconsistent with relativity theory. Similarly, the flux density of 
radiation energy transforms according to the relation S' = S (1 — 2 v ■ e/c) and not according to 
S' = S (1 — v ■ e/c) which corresponds to "fundamentally classical" result of Burns et al. (1979, p. 
5). Similarly, the fundamental condition is p' 'outgoing = (1 - Q'pr I Q'ext) p' incoming (see Eq. 16) 
and not the "fundamentally classical" condition p' outgoing = (1 — Q'pr) P' 'incoming corresponding 
to ideas of Poynting (1903), Robertson (1937), Wyatt and Whipple (1950), Burns et al. (1979, pp. 
10-11) and others. Similarly, Eqs. (7) and (11) hold. Similarly, if only thermal emission is considered, 
then the decrease of particle's mass is dm/dr = — E' Q /c 2 , where E' Q is the outgoing energy per unit 
time (see Eqs. 40-41 in Klacka 2008b). Similarly, light does not obey the law of addition of velocities, 
according to Einstein's second postulate ["Principle of the invariance (constancy) of the speed of light: 
The speed of light in a vacuum has the same value in each inertial reference frame, irrespective of 
the velocities of the light source or the light receiver." (Freund 2008, p. 9)] Simply, the processes 
participating in the P-R effect are correctly explained using ideas of relativity theory and not ideas of 
classical nonrelativistic physics. 

6.1. Burns et al. (1979) 

Since the paper by Burns et al. (1979) is generally considered to be the most relevant contribution 
on the P-R effect, we will present the most important comments on the paper. Some information can 
be found also in papers by Klacka (1992a - Sec. 2.5, 2008a, 2008b). 

1. 

Burns et al. (1979, p. 5) write: 

"The force on a perfectly absorbing particle due to solar (photon) radiation can be viewed as 
composed of two parts: (a) a radiation pressure term, which is due to the initial interception by the 
particle of the incident momentum in the beam, and (b) a mass-loading drag, which is due to the 
effective rate of mass loss from the moving particle as it continuously reradiates the incident energy." 



12 

Physics: 



Klacka et al.: Poynting- Robertson effect 



We have Q' ext = 2, Q' pr = 1 for the case (see Sees. 7.3.3 and 7.3.4 in Klacka 2008b). As for the 
incoming radiation, we can write, on the basis of Eq. (1) (see also Eq. 75 in Klacka 2008b), 
dvtJ dT = ( 2 w 2 SA'/c) W , 
(dm/dr) in = 2 w 2 SA'/c 2 , 
or, to the first order in v/c, 
dp in /dt = (2 SA'/c) (1 - v ■ e/c) e . 

The outgoing radiation effects on the particle by the following force (Eq. 76 in Klacka 2008b): 
(dp^ ut /dr) particle = - (w 2 SA'/c) {b» + u"/c) , 
(dm/dT) out = - 2 w 2 SA'/c 2 , 
or, to the first order in v/c, 

{dpout/dt) pa rticle = - (SA'/c) [ (1 - v ■ e/c) e + v/c } . 

If (2 SA'/c) (1 — v ■ e/c) e would be the "radiation pressure term", then — (SA'/c) [ (1 — v ■ e/c) e 
+ v/c] should be the "P-R drag" ("mass loading drag")! The readers should be utterly confused by 
the "physics" (see also Sec. 6.5; probably, the authors Burns et al. had not in mind these equations, 
when they formulated the above presented text). 

2. 

Burns et al. (1979, p. 5) write: 
"5' = S (1 - r/c) , (1-B) 

... The momentum removed per second from the incident beam as seen by the particle is then (S'A/c) 
S." 

Physics: 

The momentum removed per second from the incident beam as seen by the particle is (S'A/c) S , 
where S' = S (1 - 2 r/c) , 

if the accuracy to the first order in r/c is used (see also Sec. 2.5 in Klacka 1992a, or, Eqs. 20-25 in 
Klacka 2004, or, Klacka 2008b). 

3. 

Burns et al. (1979, p. 5) write about the perfectly absorbing particle: 
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"The absorbed energy flux S'A is continuously reradiated from the particle. Since the reradiation is 
nearly isotropic (small particles being effectively isothermal), there is no net force exerted thereby 
on the particle in its own frame. However, the reradiation is equivalent to a mass loss rate of S'A/c 2 
from the moving particle - which has velocity v as seen from the inertial frame of the Sun - and, 
as measured in the solar frame, this gives rise to a momentum flux from the particle — (S'A/c 2 ) v, 
schematically shown in Fig. 2." 

Physics: 

The force exerted on the particle, in its own frame, due to the outgoing radiation, is: 
(dpZ ut /dT) particle = - (w 2 SA'Q' ext /c) {<?' u"/c + (1 - Q>) &"} , 

Q = Q'pr/Q'ext J 

or, putting v = and Q' pr = 1, Q' ext = 2, 
(dp'out/dr) particle = - (S'A'/c) e' , 

according to Eq. (1). Thus, the statement of Burns et al. (1979, p. 5) that "there is no net force 
exerted thereby on the particle in its own frame" is not true. 

The following statement of Burns et al. (1979, p. 5) is: "However, the reradiation is equivalent to 
a mass loss rate of S'A/c 2 from the moving particle". We have, in reality, according to Eq. (11), 
{dm/d T ) out = -2 w 2 SA'/c 2 . 

since Q' ext — 2 for the case discussed by the authors. 

Again, the authors write: "this gives rise to a momentum flux from the particle — (S'A/c 2 ) t>, 
schematically shown in Fig. 2." In reality, the force exerted on the particle due to the outgoing 
radiation is, according to Eq. (1): 

(dp^ ut /dT) particle = - (w 2 SA'Q' ext /c) {Q' u^/c + (1 - Q') , 

Q' = Q'pr/Q' ext > 

or, for the case treated by the authors, 
(dp^ ut /dT) particle = - (w 2 SA'/c) (u^/c + &") . 
Approximation to terms of order v/c yields: 
(dp ou t/dt)p arUcle = - (w 2 SA'/c) (v/c + e/w) . 

4. 

On the basis of the correct results summarized and discussed in the previous three points (1., 2. and 
3.) the reader immediately see that also Fig. 2 is physically incorrect. The outgoing radiation in the 



14 Klacka et al.: Poynting- Robertson effect 

particle's frame is not isotropic. Also the results S' = (1 — 2 r/c) S are not taken into account. Also 
the result 

{dp ou t/dt) particle = - (w 2 S A' / c) (v/c + e/w) 
is not considered. 

5. 

Burns et al. (1979, p. 6) write: 

" The net force on the particle is then the sum of the forces due to the impulse exerted by the incident 
beam and the momentum density loss from the particle or, for a particle of mass m, 
mv = (S'A/c) S - (S'A/c 2 ) v ... (2)" 

Physics: 

The equation does not corresponds to correct physics. 

At first, since 5" is the integrated flux density measured by the moving particle (Burns et al. 1979, 
p. 5), the quantity (S'A/c) cannot be used as the relevant quantity for the equation of motion of 
the particle in the rest frame of reference of the source of light (Sun) . The second Newton's law does 
not admit such kind of mixing between two different reference frames (simultaneous considering of 
momenta measured in two different reference frames). 

Secondly, the statement of Burns et al. (1979 - Eq. 2) that 
dp in /dt = (S'A/c) S 

(dp ou t/dt) par ticle = ~ (S'A/c 2 ) V 

is incorrect. The forces acting on the particle are given by Eqs. (1), or, to the first order in v/c: 
dpi n /dt = (2 w 2 SA'/c) e/w , 

(dp ou t/dt) particle = - (w 2 S A' / c) (v / c + e/w) , 
w = 1 — v ■ e/c . 

6. 

Burns et al. (1979, p. 8) write: 
" ... we have that 

Frad = m V = Q pr M' (c - v) , 

where M' is the "mass" of the photons which strike the particle or (S'A/c 2 ). In the above form, the 
radiation force can be viewed as a drag force caused by the relative velocity c — v, of the particle 
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through a beam of photons. Finally, to terms of order v/c, we can write the net force on the particle 
as 

m v = {SA/c) Q pr [ (1 - r/c) S - v/c] . (5) 

For a totally absorbing particle this reduces to (2) since then Q pr = 1. Equation (5) is identical to 
Robertson's expression except for being more general by the inclusion of the important factor Q pr . We 
note, however, from the above derivation, that radiation forces - including the Poynting-Robertson 
drag" are fundamentally classical forces and arc not produced by relativistic effects as commonly 
believed." 

Physics: 

Our Eq. (8) yields that 
M' = (dm/dr) in = w 2 SA'Q' ext /c 2 . 

Even if we take into account that S' = w 2 S (although we have already discussed that Burns et al. 
use incorrect result given by their Eq. 1), the correct result for M' differs from the result presented 
by the authors. Substituting the correct result for M' into the equation presented above Eq. (5) of 
Burns et al, the authors would not obtain Eq. (5). 

We must admit that we have received the correct result for M' on the basis of relativity theory 
(Klacka 2008a, 2008b) and not on the basis of fundamentally classical physics, as was the access of 
the Burns et al.. 

There is another problem with the access of the authors. They use "the relative velocity c — v, of 
the particle through a beam of photons" . There is a problem. We know that we can write S' — n' h 
v' c and, similarly, S = n h v c, in the case of relativity theory (Klacka 2008b - Eqs. 20-23). Using 
transformations between n and n', v and v 1 ', we obtain the correct relation S' = w 2 S. One can say 
that the transformations between n and n' , v and v' can be obtained also on the basis of fundamentally 
classical physics (see Sec. 2.4 in Klacka 1992a and Doppler effect). And the result is consistent with 
the result known from electromagnctism (see Sec 2.5 in Klacka 1992a). But the access still assumes 
that c' = c, so there is no fundamentally classical physics. How is it possible? The explanation is 
evident: the result of electromagnetism is consistent with Maxwell's equations and they are consistent 
with relativity (Maxwell's equations motivated the creation of relativity) . If one would like to use also 
the transformation suggested by the authors, "the relative velocity c — v" , one would obtain, on the 
basis of fundamentally classical physics, that S' = w 3 S. And this result would be different from the 
correct S' = w 2 S and also from the result used by Burns et al. (1979 - Eq. 1) S' = w S. 
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It seems to us that the statement " radiation forces are fundamentally classical forces and are not 
produced by relativistic effects as commonly believed" (Burns et al. 1979, p. 8) is irrelevant. The 
complete correct results were obtained on the basis of relativistic approach (Klacka 2008a, 2008b) 
and the results are, at least in some partial steps, different from those presented before the two 
papers (Klacka 2008a, 2008b). 

7. 

As for the relativistic approach presented by Burns et al. (1979, Sec. IV), the authors consider 

Po = (1 - Qpr) E t /c ... 

(see Eqs. 10, 11, 14 in Burns et al. 1979). Although the motivation for the value 
Ei = (l- r/c) SA 

is not clear (see also Klacka 1992a - Sec. 2.5), the important result is that the usage of p = (1 
— Qpr) •■• would produce inconsistencies between relativity theory and Mie's solution of Maxwell's 
equations (Klacka 2008a, 2008b). The approach of Burns et al. (1979) leads to the inequality < 
Q pr < 2, which is not consistent with the Mie's solution of Maxwell's equations. Thus, even the 
relativistic approach by Burns et al. (1979, Sec. IV) is not physically correct. It is not consistent with 
relativity theory. The correct equation is given by Eq. (16). Moreover, the correct equation for Ei is 
Ei = wSA'Q' ext 

(Klacka 2008b - Eqs. 11, 19, and 23) and this is different from the result presented by Burns et al. 
(1979 - below Eq. 16). 

8. 

We have discussed the result for (dm/d,T) in in Sec. 2.1. The approach of Burns et al. (1979) 
would yield {dm/dr) in — wSA' /c 2 (but the authors did not treat the problem in a relativistically 
covariant form, so they could not present the result for (dm/dr) in ). Thus, the velocity dependence 
would be given by w 1 and not by w 2 . Moreover, optical properties of the particle would be repre- 
sented only by the particle geometrical cross-section A' and not by the extinction cross-section A'Q' ext . 

9. 

Burns et al. (1979, Sec. V) discuss "solar wind corpuscular forces" . We want to say that the right-hand 
side of their Eq. (17*) must be multiplied by a factor x' depending on material properties of dust 
particle (1 < x' < 3, approximately, see Eq. 29 in Klacka and Saniga 1993). The force for the solar 
wind (if only radial velocity of the solar wind is considered) is: 
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F sw = F sw [ (1 - v ■ e/v sw ) e - x' v/v sw ] , 

where F sw is the force on the particle and v sw is the heliocentric solar-wind speed. 
6.2. Shive and Weber (1982), Dohnanyi (1978) 

Shive and Weber (1982, p. 221) present: "The particle absorbs light and heat from the Sun and re- 
radiates photons uniformly in all directions. However, as seen by an observer in a frame of reference 
stationary with respect to the Sun, the photons emitted by the particle in the forward direction of its 
orbital motion will be Doppler-shifted towards higher frequencies than those emitted in the backward 
direction. The energy density ahead of the particle will accordingly be greater than that behind it, 
and the particle experiences a net tangential radiation pressure which acts as a drag on its orbital 
motion. So it spirals inward and is eventually engulfed by the Sun. Most of the interplanetary dust 
originally present at the creation of the solar system has long ago been swept out by this process." 

Similar explanation can be found in Dohnanyi (1978, pp. 562-563): "In the system of coordinates 
in which the particle is at rest, the Sun's rays (and the force of radiation pressure as well) will not 
appear to come from the centre of the Sun, but at an angle and in such a direction as to oppose 
the particle's transverse component of motion. Analysis of this effect, known as the aberration of 
light ... The opposing transverse force is then [equal to the P-R effect]. An alternative explanation of 
this effect is obtained when we view the particle's motion from the system of coordinates in which 
the Sun is at rest. The direction of the Sun's rays and that of the force of radiation pressure on 
the particle are now exactly in the radial direction. Because of the particle transverse motion, the 
radiation reflected, diffracted, or emitted by the particle is shifted in its wavelengths asymmetrically: 
radiation leaving the particle in the direction of motion is blue-shifted and in the other direction it is 
red-shifted. The net result is a drag force of the same magnitude as equation [for the P-R effect]." 

Physics of the phenomenon: 

1. 

If the case treated by Poynting (1903) and Robertson (1937) is considered, then Q' pr /Q' ext = 1/2 and 
Eq. (12), or Eq. (20), yields that the outgoing radiation causes acceleration of the particle and not its 
deceleration (dv/dt oc + v/c). This is different from the statement done by Shive and Weber (1982) 
and Dohnanyi (1982). If Q' pr / Q' ext — 1, then acceleration of the particle equals to zero in both 
reference frames, i. e., in the frame of reference of the particle and in the frame of reference stationary 
with respect to the Sun (see Eqs. 12 or 20). The important phenomenon is that mass of the particle 
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decreases in the process of reradiation and this leads to an increase of particle's acceleration: dp/dt 
= dv/dt = (1/m) [ — (dm/dt) v + ...] and the value dm/dt is given by Eq. (11) in the case of the 
P-R effect. 

Maybe, there exists more simple illustration of the correct physics. If the particle is not irradiated 
and only outgoing radiation due to the thermal emission exists, then Eqs. (40)-(41) in Klacka (2008b) 
hold. As for the spherical particle, one obtains 



acting against the motion of the particle exists, but the second of Eqs. (29) states that no acceleration 
of the particle exists and this holds in all reference frames. The particle is decelerated if the following 
asymmetry exists, in the proper frame of the particle: the total energy per unit time emitted in 
the forward direction of particle's motion (with respect to the surroundings) is greater than the 
corresponding energy emitted in the backward direction (see Eqs. 40-41 in Klacka 2008b). 

The flux density of radiation energy is influenced not only by the Doppler effect, but also by the 
change of concentration of photons: S' = w 2 S. Moreover, the term —u^/c in Eq. (14) is generated 
by the conservation of particle's mass (see Sec. 7.4.5 in Klacka 2008b) and not by the shifts in 
wavelengths or frequencies. Similarly, the last two equations of Eqs. (29) immediately lead to the first 
of Eqs. (29), if one uses the fact that p M = m u^. 



As for the force acting on the particle due to the outgoing radiation for the P-R effect, one has to use 
Eqs. (1) and (9): 

dp^/dr = - ( w 2 SA'Q' ext /c ) [ (Q' pr /Q' ext ) vT/c + (1 - Q' pr /Q' ext ) V } . 

It is easily seen that this force is not of the type represented by Eq. (29), if Q' pr /Q' ext 1 {Q' pr /Q'ext 
= 1/2 for the case treated by Poynting, Robertson, Shive and Weber, Dohnanyi). Thus, the result 
does not correspond to the idea of Shive and Weber (1982) and Dohnanyi (1978). 




(29) 



where E' 



is thermal emission energy per unit time. The first of Eqs. (29) really shows that a force 



2. 



3. 

Let us consider the following simple thought experiment. Photons are emitted from the two opposite 
rectangular faces of a cuboid (rectangular prism). The total outgoing energy per unit time is E' Q = 2 
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S'A', where A' is geometrical cross-section of the rectangular faces. Then 
(dp^/dr) e = - 2 (S'A'/c) u"/c , 

according to Eq. (29). This result can be obtained using the idea presented by Dohnanyi (1978) 
and Shive and Weber (1982). However, also the change of concentration of photons has to be taken 
into account, not only the Doppler effect: 5" = w'j-Sf — w^Sb, where the subscripts "f", "b" denote 
abbreviations of "forward" and "backward", Wf — 7(1 — v/c), Wb = 7(1 + v/c). According to Eqs. 
(18)-(20) and (24) in Klacka (2008b), we can write 
(dp»/dr) e = - [ (w}S f A'/c) b) + (wlSbA'/c) ] , 

where bj = (1/vjf ; v/vjf), b^ = (1/wb ; — v/wb); v is unit vector in the forward direction and 
orientation of the emitted photons. As a consequence, deceleration of the cuboid equals to zero: 
du^/dr = 0. 

Comment: 

The four- vectors bj, b^ are important in the previous equation. We can easily show it as follows. The 
forward and backward directions would yield (to the first order in v/c) Sf = S' (1 + 2 v/c), Sb = S' 
(1-2 v/c). But the form dp/dt = - (S f A'/c - S b A'/c) = - 4 (S'A'/c) v/c would not correspond 
to the physical result. Surprisingly, the idea of Burns et al. (1979) Sf(B) = S' (1 + v/c), Sb(B) = 
S' (1 — v/c) - results of classical mechanics - would yield the correct final result! Perhaps, Shive 
and Weber (1982) were influenced by the statement of Burns et al. (1979). Maybe, classical access 
would be of the form dp/dt = - [ S f A'/(c + v) - S b A'/(c - v)] = - 2 (S'A'/c) v/c. However, 
understanding nature requires much more effort than simple collection and combination of some ad 
hoc formulae. 

4. Aberration of light? 

As for the heuristic explanation of the term containing —v/c in the P-R effect as a consequence of 
the aberration of light, we refer the reader to Sec. 7.4.5 in Klacka (2008b). Moreover, this can be 
easily seen from Eqs. (l)-(2), or from the form to the first order in v/c: 
F in = (w\SA'Q' ext /c) (l+v e/c) e . 

The form of the force Fi n does not contain the aberrational term 
(1 + v ■ e/c) e — v/c . 

Perhaps, the following idea can appear: the aberrational term exists in the acceleration a in = v in . 
Really, something like this is produced by the covariant equation Eq. (8): 
dv m /dt = [wlSA'Q' ext /(mc)] [(l + v e/c) e - v/c] 

and this term is even almost consistent with the P-R effect! However, there is the term Q' ext instead 
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of the term Q' pr , and, the relevant four-vector — u^/c is not aberration of light. And, for the 
special case treated by Poynting (1903), Robertson (1937) and others, Q' ext — 2, Q' pr — 1 (Klacka 
2008a, 2008b). Thus, aberration of light does not explain the term — v/c in the P-R effect. 

6.3. Mignard (1982) 

Mignard (1982, p. 349) writes: "Q [Q' pr ] is a global coefficient which expresses the particle surface 
properties related to the incoming light." 

Physics: 

The dimensionless efficiency factor for radiation pressure (integrated over stellar spectrum and 
calculated for radial direction) considers volume properties of the particle, in general. It results from 
the solution of Maxwell's equations (Mie 1908). 

6.4. Kapisinsky (1984) 

Kapisinsky (1984) discusses disturbing and dissipative (destructive and disruptive) effects acting on 
interplanetary dust particles. The author explains the action of the solar electromagnetic radiation 
(Kapisinsky 1984, disturbing forces - p. 104): 

"Effect of solar electromagnetic radiation (direct light pressure): 
This effect is closely related to the Poynting-Robertson effect. Among the non-gravitational effects, 
these two are of the greatest importance for studying the dynamics of small particles in the Solar 
System. The outdirected solar pressure always acts on the particle against the gravitational force of 
the Sun." 

Continuation on the same page is: 
"Poynting-Robertson effect: Apart from the outdirected solar pressure acting on the particle, the 
absorption of solar energy by the particle and its isotropic emission causes a small force to be 
generated along the tangent to the trajectory which decelerates the particle in the orbit. This 
tangential force decreases the kinetic energy and the orbital angular momentum of the particle. 
Consequently, the particle is forced to move to an orbit closer to the Sun." 



Physics: 
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The "effect of solar electromagnetic radiation" on spherical dust particle is equivalent to the 
Poynting-Robertson effect, not to the "direct light pressure". Direct light pressure is a velocity inde- 
pendent part of the P-R effect. Thus, it is not possible to classify "direct light pressure" and the P-R 
effect as two different disturbing effects influencing motion of the particle. 

Simlarly, we have the Poynting-Robertson effect and not three effects: "direct light pressure", 
"Poynting-Robertson drag" and the "Poynting-Robertson effect". We have only the "Poynting- 
Robertson effect" and the other two "effects" are indispensable parts of the "Poynting-Robertson 
effect" . The two effects cannot be physically separated from the P-R effect. Relativistically covariant 
equation of motion, represented by Eq. (14), describes the P-R effect and it contains both the "direct 
light pressure" and the " Poynting-Robertson drag" . 

Let us consider perfectly absorbing spherical particle, within geometrical optics approximation. 
The "isotropic emission" is, in reality, described by the value (see Eq. 16) 

Q' = Q'pr/Q'ext = 1/2 , 

and by the condition 

P' outgoing = ( V 2 ) P' incoming • 

Subsequently, the force corresponding to the incoming and outgoing radiation contains also a compo- 
nent tangent to the trajectory which decelerates the particle in the orbit. This tangential component of 
the force decreases the total (kinetic plus potential) energy of the particle. The tangential component 
of the force increases kinetic energy of the particle. 

6.5. Leinert and Grun (1990) 

Leinert and Grim (1990, p. 226) write: 

Seen from the moving dust particle, the incident solar radiation is displaced by an angle v tan /c from 
the radial direction, where v tan is the tangential component of the orbital velocity. Seen from outside, 
the particle motion causes scattering and thermal emission to have a forward component. The result 
is the same: a braking force, on the average over one orbit directed opposite to the direction of motion, 
equal in first order to v/c times the radiation pressure force for a particle in a circular orbit." 

Later on, the authors state: 
"The impact of solar wind ions also exerts an outward pressure on orbiting dust particles. It is more 
than three orders of magnitude weaker than radiation pressure. But the aberration angle Vtan/v SW} 
where v sw is the velocity of the solar wind, is much larger, making the resulting "ion impact" drag 
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comparable to the Poynting-Robertson effect." (Leinert and Griin 1990, pp. 226-227). 



Physics: 

1. 

At first, Leinert and Griin (1990, Sees. 5.4.2 and 5.4.3) discuss "radiation pressure" and the 
"Poynting-Robertson effect" as two different phenomena. However, the P-R effect is the effect of 
radiation pressure on a moving body. 

2. 

As for the forces, the mathematical forms for the P-R effect seem to be consistent with the statements 
of the authors. However, the force for the solar wind (under the assumption that only radial velocity 
of the solar wind is considered) is: 
F sw = F sw [ ( 1 - v ■ e/v sw ) e - x' v/v aw ] , 

where F sw is the force on the particle, v sw is the heliocentric solar-wind speed, and, the value of 
x' , 1 < x' < 3 (approximately), depends on material properties of the particle. This force does not 
correspond to the statement of the authors. The presence of the term x' shows that any idea about 
importance of the aberration effect is not acceptable (see also Sec. 7.4.5 in Klacka 2008b, as for light 
pressure). 

Moreover, one has to take into account that Eqs. (1) hold. Even if the authors want to take into 
account the simplest case treated by Poynting (1903) and Robertson (1937), the explanations of the 
authors will not produce the consistent results for the incoming and outgoing radiation: it is sufficient 
to consider that Q' ext = 2, Q pr — 1 (see Klacka 2008a, 2008b). The tacitly assumed values Q' ext = 1, 
Q' pr = 1, considered by Poynting (1903), Robertson (1937), Burns et al. (1979) and others, arc not 
correct. 

3. 

Let us look if the substitution "force" — ► "acceleration" would make the statements of the authors 
correct. 

Eqs. (8) and (20) (see also Eqs. 106-107 of Klacka 2008b) yield 
(dv/dt) in = a ph [ (1 - v ■ ejj/c) e R - v/c ] , 

and Eqs. (12) and (20) (see also Eqs. 109-110 of Klacka 2008b) yield 
(dv/dt) out = - a ph ( 1 - Q pr /Q' ext ) [ (1 - v ■ ejj/c) e R - v/c ] ; 
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a P h — SA'Q' ext /(mc). If Q'p r /Q' ex t < lj then "scattering and thermal emission" cause acceleration 
of the particle [the case treated by Poynting (1903), Robertson (1937), and others, corresponds to 
Q'pr/Qext = 1/2], (dv/dt) ou t oc + v / c. This result differs from the statement of Leinert and Griin. 
The incoming radiation causes deceleration of the particle, (dv/dt) in oc — v/c. 

6.6. Banaszkiewicz et al. (1994) 

The authors do not try to explain the P-R effect, but they use the electromagnetic Poynting-Robertson 
force F e in the form (Banaszkiewicz et al. 1994 - Eq. 10) 



Besides misprinting in the exponent of r in the denominator, the numerical factor 2 at the term r 
I c is not correct and the correct factor equals 1. 

Maybe, Banaszkiewicz et al. (1994) were aware that the flux of the incoming radiation, as seen by 
the particle, isS' = S(l — 2r/c) and not by the relation used by Burns et al. (1979, p. 5 - Eq. 
1). However, correct equation of motion is 



see also Eq. (22) in this paper. 
6.7. Gustafson (1994) 

Gustafson (1994, p. 566) writes: "Poynting-Robertson drag can be thought of as arising from an 
aberration of the sunlight as seen from the particle and a Doppler-shift induced change in momentum. 
To the first order in v/c, the radiation force acting on a spherical particle is 



The author continues: "The second term along f is due to Doppler shift. With the transverse last 
term, this velocity-dependent part of Equation (6) [Eq. 6-G in our notation] is the Poynting-Robertson 
(PR) drag. Doppler shift enters twice: once due to the rate at which energy is received and a second 



F e = [S Ar 2 I (cr 2 ) } Q pr [( 1 - 2 r / c ) r - ( r/c) v ] 



where ... 



The same relevant terms are used in Banaszkiewicz {et al. 1994 - Eq. 22). 



Physics: 



F e = [S Ar 2 I (cr 3 ) ] Q pr [( 1 - r / c ) r - ( r/c) v } , 
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time due to the reradiated and scattered radiation." 



Physics: 

As it was written in Sec. 7.4.5 in Klacka (2008b), the third term does not originate from aberration 
of light. The second term is a result of combination of the change of concentration of photons, the 
Doppler shift of incident radiation, the term 1/w present in & M (see Eqs. 2 and 14) and the conservation 
of particle's mass. The Gustafson's statement, that this term is due to the Doppler shift of the absorbed 
and re-radiated radiation, is not correct. 

6.8. Considine (1995) 

Considinc (1995, p. 2534) writes: 

"Poynting- Robertson effect. The effect upon the motion of a micrometeoroid or other very small 
particle due to the absorption and emission of radiation. The particle absorbs radiation from the sun 
only in one direction, but reradiates energy in all directions. This effect produces a drag upon the 
particle that is directly tangential to the orbit of the particle about the sun, and thus decreases the 
orbital angular momentum. Since this angular momentum varies as the square root of the orbital 
radius, so that the particle follows a spiral path directing steadily closer to the sun. Although the 
solar radiation pressure upon the particle opposes this effect, it is not sufficient to offset it completely." 

Physics: 

The effect is characterized by an equation of motion and it holds for (non-rotating) spherical body, 
irrespective of its size. In general, it is not relevant that the process of absorption occurs, since the 
effect holds also for total specular reflection (see Sec. 7.3.3 in Klacka 2008b). May be, the explanation 
presented above concentrates to the case treated by Poynting and Robertson, when perfectly absorbing 
(within geometrical optics approximation) spherical particle is considered. However, even in this special 
case the radiation is absorbed only partially and the effect of diffraction plays equally important role. In 
general, the fact that the particle "reradiates energy in all directions" does not determine the behavior 
corresponding to the Poynting-Robertson effect (see experimental results of Krauss and Wurm 2004 
and Eqs. 9-11, 38-39 in Klacka 2008b). Eq. (16) must be fulfilled for the Poynting-Robertson effect. 

The effect produces a drag upon the particle that consists both from radial term and the term 
" directly tangential to the orbit of the particle about the sun" . 
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The solar radiation pressure upon the particle does not oppose the effect, since it is an inevitable 
part of the Poynting- Robertson effect. More correctly, the Poynting- Robertson effect is the effect of 
radiation pressure force for a moving spherical body. 

6.9. Murray and Dermott (1999) 

Frequent idea about the essence of the P-R effect is that "the term — v/c in the equation of motion 
exists under the hypothesis of isotropic reemission of radiation" . However, this idea is in contradiction 
with the following statement: "P-R drag is caused by the nonuniform reemission of the sunlight that 
a particle absorbs." (Murray and Dermott 1999, p. 121). 



In order to be able to correctly understand physics of the P-R effect, we have to realize that the 
fundamental condition for the validity of the P-R effect exists: see Sec. 3 and Eq. (16). 

The essence of the P-R effect is that nonradial components of the particle's radiation pressure and 
components of the thermal force are zero (in particle's frame of reference). Eq. (16) states that the 
P-R effect holds if the total momentum (per unit time) of the outgoing radiation, integrated over the 
whole space angle, is in the direction of the incident radiation, in the particle's frame of reference. 

If the spherical particle does not absorb, e. g., in the case of specular reflection, no thermal emission 
exists and particle orbital evolution corresponds to the P-R effect - Eq. (16) is fulfilled. This result 
also illustrates that the formulation " isotropic reemission of radiation is essential for the P-R effect" is 
not correct. If at least part of the incident radiation is absorbed, then, according to Eq. (16), thermal 
emission force equals zero and reemission of the absorbed radiation is isotropic. Thus, the statement 
"P-R drag is caused by the nonuniform reemission of the sunlight that a particle absorbs" is not 
correct (in the reference frame of the particle). 

In general, the P-R effect admits both isotropic condition for the outgoing radiation 

P outgoing > 



Physics: 



Q — Qpr/Qext ~ 1 ' 



(30) 



and anisotropic 




(31) 
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also. The form of anisotropy is exactly given by Eq. (31). 
Conclusion: 

Any statement on thermal emission is not sufficient for the existence of the P-R effect, in general. 
Real conditions for the validity of the P-R effect arc presented by Eq. (16). 

6.10. Srikanth (1999) 
1. 

Srikanth (1999) introduces "the two parameter-model". However, this model does not respect the 
physics followed from Maxwell's equations. Correct physics is given by Eq. (1) (Eqs. 75-76 in Klacka 
2008b). 

2. 

Srikanth (1999, p. 231) states: 

" The principal conclusion is that the motion of a dust in circumsolar orbit is governed only by solar 
radiation absorption and not by the asymmetry of reemission, even when viewed in the rest-frame of 



Motion of a body is governed by an equation of motion. As for the force generated by the incoming 
radiation, we have, on the basis of Eq. (1) (see also Eq. 75 in Klacka 2008b) 



or, to the first order in v/c, 

dp in /dt = ( w 2 SA'Q' ext /c )(l+ve/c)e , 

where Q' ext = 2 for the case treated by Poynting (1903), Robertson (1937), or, Srikanth (1999). The 
above presented equations do not generate the P-R effect. Thus, the outgoing radiation plays equally 
important role in the P-R effect: 



(d Pout /dt) partlc i e = ( SA'Q' ext /c) {Q'v/c + (1 - Q'){1 - v-e/c) e} , 

Q' = Q' P r/Q ext 

(see Eqs. 1, 2, 19). The accelerations are given by Eqs. (8), (12), (20) (see also Eqs. 106-111 in Klacka 



the Sun. 



Physics: 



• 2 SA'Q' ext /c ) b» 




2008b). 
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3. 

Srikhant (1999, p. 233) discusses the origin of the term 1 — v ■ e/c = 1 — f/c, in the P-R effect. He 
states that the term " is a flux- modification factor and not red-shift factor" . 



Physics: 



Flux-modification factor is ( 1 — 2 f / c ), as it follows from the relation 
S' = w 2 S, w = 1 — v ■ e / c . 

This is discussed in Klacka (1992a - Sec. 2.5; 1993b; 1994; 2000a; 2004 - Sec. 3.1, Eqs. 20-23; 2008b - 
Sec. 3.1, Eqs. 20-23). The change of concentration of photons and red-shift (Doppler effect) are equally 
important. The origin of the term 1 — v ■ e / c is evident from Eqs. (1) and (2) (see also Sec. 10 in 
Klacka 2008b). 



6.11. Wool f son (2000) 

Woolfson (2000, p. 401) considers perfectly absorbing spherical particle (within geometrical optics 
approximation) of radius a and density p at a distance R from the Sun. This particle absorbs energy. 
Woolfson writes: "The total energy absorbed per unit time is 
P = Lq a 2 /(4R 2 ) . (V.l-W) 

The energy comes to the particle radially but then is re-emitted by the particle moving in its orbit 
at speed v. The radiated energy thus possesses momentum which must be taken from the particle. 
Considering the mass equivalent of the radiated energy the rate of change of momentum of the 
particle, or the tangential force on it, is 
F = L Q a 2 v/(4R 2 c 2 ) . (V.2-W) 

This force exerts a torque that changes the angular momentum of the particle at the rate (the sign 
in front of Lq is corrected) 

dh/dt = - FR = - L Q a 2 v/(4Rc 2 ) . (V.3-W) 



If the mass of the particle is m then its angular momentum is ft = m ^/GMqR so that 
dh/dt = (1/2) m ^GMq/R dR/dt = (1/2) m v dR/dt . (V.4-W) 
From (V.3-W) and (V.4-W) and expressing m in terms of a and p we find 
dR/dt = - 3LQ/(8nRac 2 p) . (V.5-W) " 
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Woolfson's argument that re-radiated energy possesses momentum which must be taken from the 
particle is irrelevant because the incident (absorbed) energy also possesses momentum and it is given 
to the particle. The author, maybe, tacitly assumes that the radially incoming radiation does not 
change the particle's orbital angular momentum and he does not take it into account. 

However, according to Woolfson, the decrease of the particle's angular momentum due to the re- 
emitted radiation would be responsible for spiralling of the particle onto the Sun. This idea is wrong. 
The correct rate of change of momentum of the particle, due to the outgoing radiation, is given by 
Eq. (19) and not by Eq. (V.2-W). 

Let us look at the rate of change of the particle angular momentum due to the incoming and the 
outgoing radiation. For the outgoing radiation we get to the first order in v/c 

dh\ d , . / dm\ ( dv 

— (mr x v) = — — r x v + mr x — 



dt } „„,„ dt V dt J out \ dt 

'hi 



SA'Q' ext ...... , SA'Q 



SA'Q 



e Z±rxv + J 2 ext (1 - Q') rxv = 

- Q> rxv , (32) 

cr 

where we used Eqs. (11) and (20). Similarly, for the incoming radiation we have (see Eqs. 7 and 20) 

dh\ ( dm\ ( dv 

— = — — r x v + mr x — 

dtj m \dt) in \dt 

= SA 'Q'^ rxv - SA '&^ rxv = . (33) 
c z c z 

We see that the outgoing radiation really decreases the angular momentum of the particle. But 
acceleration is the quantity determining the particle's motion. Eq. (20) yields that (dv/dt) ou t °c 
v, for Q' = 1/2 corresponding to perfectly absorbing spherical particle within geometrical optics 
approximation. Thus, the particle would accelerate and its heliocentric distance would increase if 
the outgoing radiation would be the only relevant process. Deceleration of the particle is caused by 
the incident radiation, see Eq. (20). We can say, from this viewpoint, that the incident radiation is 
significant in spiralling of the particle into the Sun. Eq. (32) shows that the decrease of the particle 
angular momentum, in the case of the outgoing radiation, is caused only by the decrease of the 
particle's mass (re-emission of radiation). Thus, Woolfson's calculation (Eqs. V.3-W, V.4-W, V.5-W) 
is not correct, since he assumes that the change of angular momentum is connected only with the 
change of orbital radius. On the other side, the incident radiation does not change particle's angular 
momentum (see Eq. 33). 
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Formulae to the first order in v/c will be considered. 

The relation between Eqs. (V.l-W) and (V.2-W) is based on the relativistic relation 

P = f v ■ ( 34 ) 

This relation yields 
dp E E dv 

i = ^ v + ^dt- < 35 > 

Woolfson access corresponds to E in — E out and he uses —E out instead of E. Woolfson neglects 
the second term on the right-hand side of Eq. (35). He obtains 

p w = -(Eout/c 2 ) V . 



If we would use the correct result represented by Eq. (19), then 
A' O' 

Pw = ^ « ■ (36) 

This equation should correspond to Eq. (V.2-W). 

However, one can immediately see that even the result in Eq. (36) does not correspond to the 
result presented in Eq. (V.2-W). The difference can be easily explained: Woolfson's considerations are 
consistent with those of Poynting (1903), Robertson (1937), Wyatt and Whipple (1950), Burns et al. 
(1979) and others - all of them state that Q' ext = 1, but the correct result is Q' ext = 2. 

Moreover, Eq. (36) is not correct, in reality. It differes from Eq. (21). The reason is that the second 
term in Eq. (35) cannot be neglected! The first term in Eq. (35) is negligible, in reality! We will show 
it, now. 

We have, on the basis of Eqs. (18)-(19), 



dE _ dE rn _ dE out _ , = t ^. v ■ e 
dt ~ dt dt ~ b A Q ext 



S A' Q' ext Q> . (37) 



Using E = m c 2 , Eqs. (35) and (37) yield 
dp dv 

—j— — TTl —— . 38) 

dt dt K ' 

Using v — (dv/dt)i n + (dv/dt) out , Eqs. (20) yield the result consistent with Eq. (21). But Eq. (21) 
differs from Eq. (36)! 

Surprisingly, vector products r x (dp/dt) and r x p w , using Eqs. (21) and (36), yield the same 
results if Q' ext — 1 (incorrect) is assumed! 
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Let us continue. Woolfson writes h — to \JG M R. If we consider total energy E, as we have 
just discussed, then to is a constant. But, at first, the form of the relation for angular momentum 
h, corresponding to Woolfson's idea, is to y/G M & (1 - 0) R (see Sees. 6.1 and 6.2 in Klacka 2004). 
Secondly, the relation for h assumes a circular orbit. But, in this case, Eq. (37) yields dE/dt — and 
the Woolfson's access yields p = and, thus, dR/dt = 0. But it is O.K., since R does not change 
when a particle moves in a circular orbit. Of course, this is a trivial result and one does not need the 
complicated "physics" presented by Woolfson. 

Eq. (21), or Eq. (22), is the correct equation of motion (gravity of the Sun is also added, this is 



tacitly assumed) . The size of the angular momentum h is m \JG M (1 — (3) pp, where pp is semi-latus 
rectum. Eqs. (21)-(22) yield 

dt to c A c R z 

h = TO yjGMr s (1-/3) pp . (39) 

Averaging over one orbital period (see Sec. 6.1 in Klacka 2004) yields 

*--,_£*S_*, (40, 

where ap is semimajor axis and ep is eccentricity of the elliptical orbit. We can use pp = ap (1 — e^) 
and Eq. (40) yields 

*S_- 2fl _£"iL_ w . (41) 

c a^y 1 

Eq. (41) is correct equation for secular evolution of semi-latus rectum. If can be neglected in 
comparison with 1, then (but see also Klacka and Kaufmannova 1992, 1993, Breiter and Jackson 
1998, Klacka 2001a) pp equals approximately ap and Eq. (41) yields 



ap = \Ja} 4/3^® (t t ) , (42) 
if ap o is the value of semimajor axis at the time to; this results is consistent with the Eq. (V.5-W). 



6.12. McDonnell et al. (2001) 

McDonnell et al. (2001, p. 164) write that "small particulates are subject to forces which significantly 
alter or dominate their dynamics" and also the P-R effect is included: 

"the Poynting- Robertson effect (non-radial radiation pressure due to aberration of sunlight) resulting 
in a gradual decrease in eccentricity and semimajor axis" . 



Physics: 
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If we would consider arbitrarily shaped dust particle, then also non-radial components of radiation 
pressure may exist. Spherical particle exhibits only radial component of the radiation pressure. This is 
exact physics and it is defined in the proper frame of reference of the particle. In the reference frame 
of the Sun, there is also a non-radial term coming from —v / c. But this term does not come from the 
aberration of light, see Sec. 7.4.5 in Klacka (2008b), see also Sec. 6.2, and it is a consequence of the 
constancy of particle's mass. Moreover, other physical effects are important in the P-R effect (Dopplcr 
effect, change of concentration of the incident photons) and the P-R effect cannot be reduced only to 
the term — v / c. This term would not yield the correct result for the secular decrease of eccentricity 
and semimajor axis. One has to bear in mind that the P-R effect is the effect of radiation pressure 
on moving spherical particles. Moreover, if the particles are small, then they can be expelled from the 
Solar System. 

6.13. Griin et al. (2001a) 

Grim et al. (2001a, p. 335) write: "Besides the direct effect of radiation pressure on trajectories of 
small dust grains there is also a more subtle effect: the Poynting-Robertson effect. This is caused by 
radiation pressure force, that is not perfectly radial on moving dust particle but has small component 
opposite to the particle motion. The strength is order of v t /c (v t is the tangential particle velocity) 
and it leads to loss of angular momentum and orbital energy of orbiting particle. The effect is 
strongest when the particle speed is highest, i.e. close to the Sun at its perihelion. Thefore particle 
orbits slowly circularized while they spiral toward the Sun." 

Physics: 

The Poynting-Robertson effect is the effect of radiation pressure on moving spherical particles. 
The "direct effect of radiation pressure" is an indispensable part of the P-R effect. There are not 
two different effects, "the direct effect of radiation pressure" and "a more subtle effect: the Poynting- 
Robertson effect" , as the authors state. 

The P-R effect is caused by radiation pressure force which is perfectly radial. And this holds in 
the proper frame of reference of the particle. Yes, the effect contains also a velocity component v t /c, 
which leads to the loss of angular momentum and orbital energy of the orbiting particle. But v t is 
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the transversal and not the tangential particle velocity. The velocity vector is always tangential to the 
orbit of the particle. The velocity vector can be decomposed into radial and transversal components. 



6.14. Dermott et al. (2001) 

Dermott et al. (2001, pp. 574-575) write: "The component of the radiation force tangential to a 
particle's orbit is called the P-R drag force. This force is also proportional to (3. It results in an evolu- 
tionary decrease in both the the semi-major axis and the osculating eccentricity of the particle's orbit." 

Physics: 

According to the conventional definition (see Sec. 5), the P-R drag is the part of equation of 
motion which contains velocity terms. So, not only the term containing — v/c, which is tangential to 
the particle's orbit. However, relativistic covariant formulation does not admit to decompose the P-R 
effect into various parts. 

The osculating eccentricity does not exhibit only decrease: de/dt can be also positive (see Eqs. 55 
and 94, 96 in Klacka 2004), but secular value of de/dt is negative (if analytical time averaging over 
one period is admitted). An interesting result is presented in Figs. 5-6 in Klacka et al. (2007). 

6.15. Grun et al. (2001b) 

Grim et al. (2001b, p. 791 - glossary) explain: " Poynting- Robertson drag: Drag on a moving particle 
that is due to the asymmetrical momentum distribution of scattered and absorbed light; the drag 
primarily affects small grains, causing them to spiral toward the central object that they orbit, 
whether the Sun or planet." 

Physics: 

There is a preferred axis. It is given by the unit vector of the incoming radiation Pi ncoming . Total 
outgoing radiation preserves this symmetry, since the condition given by Eqs. (16) is fulfilled: 

Poutgoing = (X ~ Qpr/Qext) P incoming ■ 

As a consequence, the equation of motion of the spherical particle 
dp'/dr = (Q' pr /Q' ext ) P incoming 
also fulfills the symmetry. 
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The drag — v/c is generated by the first condition in Eqs. (16): 
E' = E' 

outgoing incoming 

The condition corresponds to the conservation of particle's mass. 

Motion of spherical totally reflecting (within geometrical optics approximation) body is also de- 
scribed by the P-R effect. Thus, absorption of the incoming light is not an important condition for 
the P-R effect. 

6.16. Matzner (2001) 

In the book of Matzner (2001, p. 524) we can find the following explanation of the Poynting- Robertson 
effect. It is "a drag force arising on particles orbiting the Sun, because solar radiation striking the 
leading surface is blueshifted compared to that striking the following surface. Thus, the particles 
receive a component of momentum from the radiation pressure which is opposite to the direction of 
motion." 

Physics: 

If we consider a homogeneous flow of radiation from star at particle's location (it is a good ap- 
proximation, in practice), then this radiation impacts on the particle in one direction. There is no 
difference between Doppler shifts of photons striking different parts of particle's surface. The origin 
of drag force component is in the conservation of particle's mass. The role of the Doppler effect is in 
the relation S' w 2 S, where S is the flux of radiation energy (see Eqs. 1 and 2). 

6.17. Murdin (2001) 

Murdin (2001) presents that Poynting- Robertson effect is "a net force acting on small particles which 
causes their orbits to decay, also known as Poynting- Robertson drag" . His explanation of this effect 
is following: 

"When photons of light strike a dust particle in orbit around the Sun, they impinge not directly 
'side-on' but slightly on its leading side (the 'forward' side in terms of its orbital motion). Their 
energy is absorbed, and subsequently re-radiated, but it is re-radiated isotropically - in all directions. 
The transfer of energy from photon to particle can be thought of as the application of a force to 
the particle. There is a component of this force that acts in the opposite direction to the tangential 
component of the particle's orbital motion which is not compensated for when the absorbed energy 
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is re-radiated. As a result the particle loses kinetic energy, which reduces its velocity, so it spirals 
inward, toward the Sun. ... The effect can be counteracted by radiation pressure, more so with 
decreasing particle size below about 1 jum." 

Physics: 

The first mistake is that Murdin distinguishes radiation pressure and P-R drag which denotes as 
the P-R effect. In reality, only the P-R effect exists. The P-R effect is the effect of radiation pressure 
on a moving spherically symmetric body (see also Sec. 6. 2). 

It is evident that Murdin considers perfectly absorbing spherical particle. He writes that the 
incident radiation causes force with component tangential to the particle's orbit due to the aberration. 
This is not true in the stationary reference frame connected with the Sun. As we can see from Eqs. 
(18), the force originating from incident radiation has radial direction. Mentioned tangential force 
originates just from the outgoing radiation. Moreover, this tangential force decelerates the particle 
in the orbit and decreases the total (kinetic plus potential) energy of the particle. But increases the 
kinetic energy. 

6.18. Bottke et al. (2002a), Lauretta and McSween (2006) 

Bottke et al. (2002a, p. 769), Lauretta and McSween (2006, p. 914) write: " Poynting- Robertson effect 
- an effect of radiation on a small particle orbiting the Sun that causes it to spiral slowly toward 
the Sun. It occurs because the orbiting particle absorbs energy and momentum streaming radially 
outward from the Sun, but reradiates energy isotropically in its own frame of reference." 

Physics: 

The P-R effect holds for any spherically symmetric dust particle, not only for "perfectly absorb- 
ing" particle which "isotropically reemits/reradiates the absorbed radiation, in the reference frame of 
the particle". E. g., also for totally reflecting sphere (within the geometrical optics approximation), 
although no absorption of the incoming radiation occurs. Moreover, even within the geometrical optics 
approximation, one must be careful in stating that the radiation is perfectly absorbed and isotropically 
reradiated in the particle's own frame of reference. Although these ideas come back to Poynting (1903), 
and Robertson (1937), and although they are conventionally used (e. g., Wyatt and Whipple 1950, 
Burns et al. 1979), they do not contain the whole physics. In reality, the particle accepts momentum 
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per unit time given by Eq. (1), or, 
dp'Jdr = (S'A'Q' ext /c) e! , 

in the frame of reference of the particle, and, the outgoing radiation is described by Eq. (16). The 
case treated by Poynting, Robertson and others corresponds to Q' pr — 1, Q' ext = 2, in reality: 

Poutgoing = (1/2) 

incoming 

and not p' outgomg = (Poynting 1903, Robertson 1937, Wyatt and Whipple 1950, Burns et al. 1979). 

The P-R effect is the effect of radiation pressure. The pressure can cause escape of small dust 
particles from the Solar System. Only larger particles (of radius greater than 0.1 micron, approximately 
- this depends on optical properties of the particles) spiral slowly towards the Sun, due to the P-R 
effect. 

6.19. Bottke et al. (2002b), Bottke et al. (2002c) 

Bottke et al. (2002b, p. 11), Bottke et al. (2002c, p. 396) write: "Poynting-Robertson drag - a radiation 
effect that causes small objects to spiral inward as they absorb energy and momentum streaming ra- 
dially outward from the Sun and then reradiate this energy isotropically in their own reference frame." 

Physics: 

The P-R drag is not a physical effect, which can be treated as an independent phenomenon. It 
is an indispensable part of the P-R effect. It holds for any spherically symmetric dust particle. The 
"perfect absorption" is not the essence of the P-R effect, as it was discussed in the previous Sec. 6.18. 

6.20. Danby (2003) 

Danby (2003, p. 366) writes: "We shall examine the Poynting-Robertson effect. A small particle moves 
in the solar system, subject to the gravitational attraction of the Sun and to radiation pressure from 
the Sun." 

Danby continues: " ... aberration of light ... leads to a transverse drag component of the radiation 
force, proportional to the transverse component of velocity, divided by the speed of light." And then: 
"Next, consider radial motion away from the Sun with speed f. Because of this, the particle will 
absorb less radiative energy in unit time than it would if it were at rest. Also, the radiation that it 
receives is diluted by the Doppler effect. So, relative to the 'unperturbed' model, we have a radial 
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drag force proportional to 2 r/r 2 ." 



Physical comments: 

1. 

It must be said that spherical particle with spherically distributed mass is considered. Moreover, 
equations presented by the author assume that optical properties of the particle do not change. It 
may be wise to be stressed, too. 

2. 

As for the explanation of the term containing — v/c in the P-R effect as a consequence of the 
aberration of light, we refer the reader to Sec. 7.4.5 in Klacka (2008b), and, also to Sec. 6.2. The 
term is not a consequence of the aberration of light. 

3. 

Danby states, that the term — 2 f/c comes from a decrease of radiative energy and Doppler effect. 
This does not correspond to physics. The access of Danby corresponds to the procedure presented in 
Burns et al. (1979) which was commented in Sec. 2.5 in Klacka (1992a), see also Sec. 6.1: this would 
require violation of the second Newton's law, if one would like to obtain the correct formula for the 
Poynting-Robertson effect. 

4. 

According to physics, the energy flux density measured by the particle generates the term — 2 
f/c (see Eqs. 113-114 in Klacka 1992a, Eq. 23 in Klacka 2004) and the term is generated by the 
decrease of concentration of photons and frequency, i.e. Doppler effect (see Eqs. 32 and 36 in Klacka 
1992a, Eqs. 7-8 in Klacka 1993b, Eqs. 17 and 22 in Klacka 2004, Eqs. 17, 20, 22 and 23 in Klacka 2008b). 

5. 

The radial drag force proportional to 2 r/r 2 comes from the transformation (if the first order in v/c 

is taken into account) 

S' = w 2 S = (1 - 2 f/c) S , 

then from the term 

b = e/w = (1 + f/c) e , 
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and, finally, from the radial part of the term 
- {v/c) rad = - r/c . 

The term 1/r 2 comes from the radiation energy flux density. 
6.21. Parker (2003), Moore (2002), Kohler and Mann (2002) 

Dictionary of Astronomy (Parker 2003, p. 99) explains the Poynting-Robcrtson effect in the following 
way: "The gradual decrease in orbital velocity of a small particle such as a micrometeorite in orbit 
about the sun due to the absorption and reemission of radiant energy by the particle." 

Similar explanation is presented in Moore (2002): 
" Poynting-Robertson effect. Non-gravitational force produced by the action of solar radiation on 
small particles in the Solar system; it causes the particles to spiral inwards towards the Sun. When 
particles in orbit around the Sun absorb energy from the Sun and rc-radiate it, they lose kinetic 
energy and their orbital radius shrinks slightly. The effect is most marked for particles of a few 
micrometers in size." 

Similarly, Kohler and Mann (2002) explain: 
"The Poynting-Robertson effect reduces the kinetic energy of the particles. The semimajor axis and 
the eccentricity of the orbits decrease so that finally the particles fall into the stars." 

Physical explanation is based on the two-body problem: 
At first, we immediately formulate the correct statement: the P-R effect reduces the total (kinetic 
plus potential) energy of the particles. 

The Keplerian orbital velocity of a body orbiting the Sun is characterized by magnitude 



where a and r are semi-major axis and heliocentric distance of the body of mass m. Time averaging 
over one period yields 



Let us consider the P-R effect and let us assume that (3 is a constant. Then, M sun + m —> M sun 
for a micrometeoroid. One possibility is to use as a central Keplerian acceleration — GM sun r /r 3 and 




(43) 




G (M sun + to) 



(44) 
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the corresponding secular value of semi-major axis is a. If the non- velocity radiation term is included 
into the central acceleration, then the central Keplerian acceleration is — GM sun (l — (3) r/r 3 and the 
corresponding secular value of semi-major axis is ap. We can write 

{V >) = GM ^~P) (45) 
ap 

for a micrometeoroid. Since semi-major axis is a decreasing function of time for the P-R effect, also 
Jt ^ > ° ' ™ 
Qualitative physical explanation: 

The above cited Parker's statement would be correct under the assumption that potential energy 

is constant. If a particle is decelerated, then its total energy is a decreasing function of time. The time 

averaged values of the total, potential and the kinetic energies (per unit mass) of the particle are 

{E) = _ GM sun (1 - g) 
2 ap 

(U) = 2{E), 

(K) = -(E). (47) 

We want to stress that Eqs. (47) are results corresponding to secular evolution. The last two equations 
of Eqs. (47) are similar to the virial theorem, where energy E is constant. Eqs. (47) immediately yield 
d GM sun (l-P) dap 

Jt {E) ~ ~~aT ' 

d_ = GM sun (1-/3) dap Q 
dt [ ' a 2 p dt 

The particle's potential energy decreases more rapidly than the total energy, and, thus, the particle's 
kinetic energy must increase in order to reduce the decrease of potential energy. The increase of (K) 
= (v 2 )/2 occurs. So, the particle's velocity is an increasing function of time. 



Quantitative physical explanation: 



Squared orbital velocity, averaged over one period, is 
, _ GM sun (l-0) e% 1-ej 

V ' ~~ cir 1 - P 2 4 / 5 ' 1 ' 
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where ag in and ep in are initial semi-major axis and eccentricity of the particle. The value of the 
eccentricity for a time t is given by the following differential equation: 



Eq. (50) shows that eccentricity is a decreasing function of time. As a consequence, (v 2 ) is an increasing 
function of time, according to Eq. (49). 

Finally, we have to stress that the P-R effect holds also for nonabsorbing particles (totally reflecting 
sphere). Although the case treated by Poynting (1903) and Robertson (1937) holds for perfectly 
absorbing nonrotating spherical particle, the condition for absorption is not relevant (and, moreover, 
Eq. 16 is valid). 

6.22. Encrenaz et al. (2004) 

Encrcnaz et al. (2004, pp. 189-190) write: 

"The particles in the rings are subject to two effects: 

- solar radiation pressure, which pushes the particles towards the outer regions of the Solar System. 
Because it is proportional to the surface area of a grain, while gravitational attraction is propor- 
tional to the volume, radiation pressure is particularly effective on small-sized particles, about one 
micrometre across; 

- the Poynting-Robcrtson effect results from the fact that a particle orbiting the Sun receives 
solar radiation proprotional to its cross-section, but radiates it away isotropically. The result is 
preferentially radiation in the forward direction (the frequency of the photons is increased because 
of the velocity of the particle) and this translates to a loss energy by the particle. This form of 
braking, which continuously decreases the eccentricity of the grain's elliptical path, tends to turn 
the latter into a circle, and finally into a spiral. It is calculated that a grain one micrometre across, 
subject to the Poynting-Robertson effect alone would fall into the Sun in a few thousand years. In 
practice, grains of micrometre size are more sensitive to the radiation pressure that tends to push 
them towards the outer edge of the Solar System. The Poynting-Robertson effect is thus particularly 




t = xT , < x < 1 



(50) 



and T is the time of spiralling of the particle into the Sun: 




(51) 
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important for particles in the centimetre size range. 



Physics: 

Radiation pressure is an indispensable part of the P-R effect. 

The incoming radiation generates the force acting on the spherical particle. The force is propor- 
tional to the extinction cross-section of the particle (see Eq. 1). If the particle is perfectly absorbing, 
within geometrical optics approximation, then p' outgoing = (1/2) p 'incoming and not p' outgoing = 
(isotropic outgoing radiation), as it was pointed out by Klacka (2008a, 2008b). 

Momentum (per unit time) of the outgoing radiation is characterized by Eq. (1), or, Eq. (19), 
where Q' — 1/2 for the case treated by Encrenaz et al. (2004): 
dp ut/dt = (SA'/c) [ v/c + (1 - v- e/c) e ] . 

A loss of energy by the particle, due to the outgoing radiation, follows from Eq. (19): 
dE/dt = -2 S A' [ 1 - (3/2) ve/c] . 

6.23. Encrenaz et al. (2004) 

Encrenaz et al. (2004, pp. 446) write: 

"The Poynting-Robertson effect must be taken into account. As explained earlier ... this effect results 
from the fact that the - radial - force produced by the photons, which propagate with the velocity 
of light c, is exerted on a body moving at a velocity v, which is not co-linear with c. The transfer of 
momentum creates a braking force, which leads to a path that spirals in towards the Sun." 

Physics: 

The Poynting-Robertson effect results from the fact that a particle moves with respect to the 
source of radiation. The final equation of motion - represented by Eqs. (14), (21), or, (22) - holds for 
any relative direction and orientation of v and e, including their co-linearity. 

6.24. Quinn (2005) 

Quinn (2005, pp. 194-195) distinguishes the radiation pressure and the Poynting-Robertson drag. 
He writes: "... as radiation is absorbed and re-radiated by the particle, the re-radiated radiation is 
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anisotropic in the rest frame of the sun. This results in a drag force opposite the direction of motion 
of the particle referred to as Poynting-Robertson Drag. The total radiation force including this effect 
to the first order in v/c is: 

Frad = (L Q Q pr A)/(4ircr 2 ) [ (l-2v r /c)r - (v e /c) 

where L© is the solar luminosity, r is the radius of the particle, and v r and vq are the radial and 
tangential component respectively of the particle velocity in the frame of the Sun. The first term in 
this expression is the radiation pressure, and the second and third terms are the Poynting-Robertson 
drag. There is a factor of 2 in the second radial term because of the combination of Doppler shift in 
absorbing radiation and Doppler shift in emission." 



Physics: 



The factor 2 in the second radial term is a result of combination of the change of concentration of 
photons, the Doppler shift of the incident radiation, the term l/w present in 6^ (see Eqs. 2 and 14) 
and the conservation of particle's mass. The Quinn's statement that the second radial term is due to 
the Doppler shift of the absorbed and re-radiated radiation, is not correct. 

Quinn's statements are in practice the same as Gustafson's ones (Gustafson 1994). Thus, the 
correct explanation of the origin of individual terms in expression for drag force is also the same. See 
also Sec. 6.7. 



6.25. Harwit (2006) 

Harwit (2006, p. 174-175; see also Harwit 1988, p. 176-177, Harwit 1973, p. 176-177) belongs to the 
authors who try to explain the Poynting-Robertson drag through a decrease of the angular momentum. 

Harwit writes: 

"Consider a grain of dust in interplanetary space. As it orbits the Sun it absorbs sunlight, and re-emits 
this energy isotropically. We can view this process from two different viewpoints. 

(a) Seen from the Sun, a grain with mass m absorbs light coming radially from the Sun and re-emits 
it isotropically in its own rest frame. A re-emitted photon carries off angular momentum proportional: 
(i) to its equivalent mass hv/c 2 , (ii) to the velocity of the grain RQ; and (iii) to the grain's distance 
from the Sun R. Considering only terms linear in V/c, and neglecting any higher terms we see that 
the grain loses orbital angular momentum L about the Sun at a rate 
dL = {hv/c 2 ) 9R 2 , (1/L) dL = {hv)/{mc 2 ) , (5 - 45 - H) 
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for each photon whose energy is absorbed and re-emitted, or isotropically scattered in grain's rest 
frame. 

(b) Seen from the grain, radiation from the Sun arrives at an aberrated angle 9' from the direction 
of motion, instead of at 9 = 270°. Hence, 

cos (9 = \cos9' + (V/c)]/[l + (V/c) cos9'} = 0, cos <9' = - V/c . (5 — 46 — H) 

Here V is 9R 1 the grain's orbital velocity, and the photon imparts an angular momentum pRcos9' = 

- (hv/c 2 )R 2 9 to the grain. 

For a grain with cross-section a g 
dL/dt = - (L sun (T g )/(4TrR 2 mc 2 ) L , (5 - 47 - H) 
where L sun is the solar luminosity. 

Either way, the grain's velocity decreases on just absorbing sunlight. From the first viewpoint, 
this happens because the grain gains mass, which it then loses on re-emission; from the second, it is 
because the grain is slowed down by the transfer of angular momentum." 

Physics: 

We have already explained in Sec. 6.11 that although the orbital angular momentum of the per- 
fectly absorbing particle decreases due to the outgoing radiation (Eqs. 32, 39), the process increases 
acceleration of the particle (see Eq. 12, Q' = 1/2). The outgoing radiation acts on the particle by 
the force containing also the term —v/c (see Eqs. 19, where Q' = 1/2). Thus, the outgoing radiation 
participates on the Poynting-Robertson drag force, but the corresponding acceleration is positive (see 
Eqs. 20)! 

Harwit's considerations are strictly based on the idea of Poynting (1903), Robertson (1937) and 
others, when the authors assume that the outgoing radiation corresponds to the isotropically re- 
cmittcd radiation. We know that this is not true. 

Moreover, even if the outgoing radiation would be equivalent to the re-emitted radiation, the 
considerations of Harwit would be incorrect. The Harwit's statement that each re-emitted photon 
carries off the same angular momentum is not true, in a stationary reference frame. One photon, 
emitted in a direction of unit vector n', carries a momentum dp' — (hv'/c)n', in the grain's proper 
reference frame. Each photon emitted in an arbitrary direction is characterized by the frequency v' . 
Thus, in the grain's reference frame, the photons are emitted isotropically and the total momentum 
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carried off by the photons per unit time is zero (F' e = 0). The momentum dp' transforms to the 
stationary reference frame as 

(52) 



hv hv | , 
dp — — n < n + 

c c 



, , . v ■ n' 7 

7-1 —5- + 1 

v z c 



dL = R x dp = — R (e x n' + e x -) , (53) 
c V c / 



due to the Doppler effect and the aberration of light. Thus, the angular momentum carried off by the 
given photon is 

h/ 
c 

if the first order in v/c is considered. 

In his second explanation (from the view of the grain), Harwit mixes terms from stationary and 
proper reference frames. The incoming radiation really hits the grain in an aberrated direction and 
the radiation acts on the grain by a force (see Eqs. 18 and 20). But particle's velocity equals zero at 
every moment, in the particle's proper frame of reference. Thus, the particle has no orbital angular 
momentum which could be decreased in the proper frame. 



6.26. de Pater and Lissauer (2006) 

de Pater and Lissauer (2006, pp. 35) state: 

"Sections ... above describe the gravitational interactions between the Sun, planets and moons. Solar 
radiation, which provides an important force for small (< 1 m) particles in the Solar System, has been 
ignored. Three effects can be distinguished: 

— Radiation pressure, which pushes particles (primarily micrometer-sized dust) outwards from the 
Sun. 

— Poynting-Robertson drag, which causes centimeter-sized particles to spiral inward towards the Sun. 

— The Yarkovski effect, which changes the orbits of meter to kilometer-sized objects due to uneven 
temperature distributions at their surfaces. 

The solar wind produces a corpuscular drag similar in form to the Poynting-Robertson drag; 
corpuscular drag is more important for submicrometer particles." 

The authors continue with a subsection "Radiation Force (micrometer-sized particles)": 

"The Sun's radiation exerts a repulsive force, F ra d, on all bodies in our Solar System. This force 
is given by: 

F rad w [ L Q A/{A-Kcr%) ] Q pr r , (2.45-P-L) 
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where A is the particle's geometrical cross-section, L Q the solar luminosity, r Q the heliocentric distance, 
c is the speed of light, and Q pr the radiation pressure coefficient. The radiation pressure coefficient 
accounts for both absorption and scattering and is equal to unity for a perfectly absorbing particle. 
Rclativistic effects produced by the Dopplcr shift between the frame of the Sun and that of the particle 
are generally small, and have been omitted from equation (2.45-P-L), but they will be considered in 
Sec. 2.7.2-P-L. The parameter is defined as the ratio between the forces due to the radiation pressure 
and the Sun's gravity: 

= | F rad I F g | = 5.7 x 1(T 5 Q pr I (pR) , (2.46-P-L) 

with the particle's radius, R, in cm and its density, p, in g cm~ 3 . Note that is independent of heliocen- 
tric distance and that the solar radiation force is only important for micrometer and submicrometer- 
sized particles. Extremely small particles are not strongly affected by radiation pressure, because Q pr 
decreases as the particle radius drops below the (visible wavelength) peak in the solar spectrum (Fig. 
2.12-P-L). The magnitude of the Sun's effective gravitational attraction is given by: 
F g ,ef f = -(1-0) GtuMq I r| . (2.47-P-L) 

It is clear that small particles with > 1 are repelled by the Sun's radiation, and thus quickly escape 
the Solar System, unless they are gravitationally bound to one of the planets. Dust released at a 
keplerian velocity from bodies on circular orbits is ejected from the Solar System if > 0.5; critical 
values of for dust released from bodies on eccentric orbits are calculated in Problem 2.32-P-L. 

The importance of solar radiation pressure can, for example, be seen in comets (Sec. I0.4.I.-P-L): 
Cometary tails always point in the anti-solar direction due to the Sun's radiation pressure. The dust 
tails are curved rather than straight as a result of the continuous ejection of dust grains from the 
comet, which itself is on an elliptical orbit around the Sun." 

de Pater and Lissauer (2006, pp. 36-37) continue with a new subsection " Poynting- Robertson Drag 
(centimeter-sized grains)": 

"A particle in orbit around the Sun absorbs solar radiation and reradiates the energy isotropically 
in its own frame. The particle thereby preferentially radiates (and loses momentum) in the forward 
direction in the inertial frame of the Sun (Fig. 2.13-P-L). This leads to a decrease in the particle's 
energy and angular momentum and causes dust in bound orbits to spiral sunward. This effect is called 
Poynting- Robertson drag. 

Let us consider a perfectly absorbing, rapidly rotating, dust grain. The flux of solar radiation 
absorbed by the grain is equal to 
[ L Q A I (4 7T r|) ) (1 - v r I c) , (2.48a-P-L) 
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where v r = v r is the radial component of the particle's velocity (i.e., the component which is parallel 
to the incident beam of light). The second term in expression (2.48a- P-L) accounts for the Dopplcr 
shift between the Sun's rest frame and that of the particle; the transverse Doppler shift is of order 
(vg/c) 2 -C 1 and will be ignored here. The absorbed flux is reradiated isotropically and can be written 
as a mass loss rate in the particle's frame of motion (using E — mc 2 ): 
[ L & A I (4 c 2 r|) ] (1 - v r I c) . (2.48b-P-L) 

As the particle moves relative to the Sun with velocity v, there is a momentum flux from the particle 
as seen in the rest frame of the Sun, since the particle emits more momentum in the forward direction 
than in the backward direction (Fig. 2.13-P-L). The momentum flux is equal to 

[ - L Q A I (4 7T c 2 r%) } (1 - v r I c) v , (2.48c-P-L) 

which can be generalized to the case in which the particle reflects and/or scatters some of the radiation 
impingent upon it via multiplication by Q pr . The net force on the particle in this more general case 
is given by: 

Frad = [ Lq Q pr A / (4 7T C r|) ] (1 - V r / c) f 

- { L Q Q pr A v I (4 vr c 2 r%) ] (1 - v r / c) v (2.49a-P-L) 
w [ L Q Q pr A I (4tt c r| ) ] [ (1 - 2 v r / c ) r - {v e / c) 9 ] . (2.49b-P-L) 
The first term in equation (2.49b-P-L) is that due to radiation pressure and the second and third 
terms (those involving the velocity of the particle) represent the Poynting- Robertson drag. 

From the above discussion, it is clear that small dust grains in the interplanetary medium disap- 
pear, with (sub-)micrometer-sized grains being blown out of the Solar System, while ccntimetcr-sizcd 
particles spiral inward towards the Sun." 

Physics: 

The authors distinguish the two effects: the radiation pressure and the Poynting-Robertson drag. 
However, physics states that the effect of radiation pressure on a moving spherical particle is equiv- 
alent to the Poynting-Robertson effect. The Poynting-Robertson effect is the relevant physical effect. 
The mentioned two " effects" , the radiation pressure and the Poynting-Robertson drag cannot be sep- 
arated from the P-R effect, as it follows from the relativistic formulation of the equation of motion. 
Similarly, the two "effects" are relevant for the same particles: it is not correct to say that one part 
of the P-R effect is relevant for the micrometer-sized dust, while the rest part is relevant for the 
centimeter-sized particles. Moreover, the P-R effect is not the dominant effect for centimeter-sized 
interplanetary particles - collisions among particles, impact erosion are important. Moreover, the P-R 
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effect is dominant in the inner part of the Solar System (say, for heliocentric distances less than 10-20 
AU - solar wind bombardment, effect of the interstellar hydrogen are relevant for larger heliocentric 
distances) for micrometer-sized (and tens of micrometer-sized) dust. Small particles can be blown out 
from the Solar System, due to the radiation pressure. 

As for the solar wind effect ("corpuscular drag"), it is important for the same sizes of the particles 
as the P-R effect. As for submicromcter-sized particles, the Lorentz force, due to the presence of 
interplanetary magnetic field, is relevant. Small particles can be blown out from the Solar System, 
due to the (electromagnetic) radiation pressure. 

As for the Sec. 2.7.2-P-L: "A particle in orbit around the Sun absorbs solar radiation and 
reradiates the energy isotropically in its own frame. The particle thereby preferentially radiates (and 
loses momentum) in the forward direction in the inertial frame of the Sun (Fig. 2.13-P-L). This leads 
to a decrease in the particle's energy and angular momentum and causes dust in bound orbits to 
spiral sunward. This effect is called Poynting- Robertson drag." 

As the authors state, the derivations presented in the textbook is based on the derivation 
presented by Burns et al. (1979). Unfortunately, the derivation is physically incorrect (Klacka 1992a 
- see mainly Sec. 2.5). We will present the relevant comments, now. 

1. 

"Let us consider a perfectly absorbing, rapidly rotating, dust grain. The flux of solar radiation 

absorbed by the grain is equal to 

[ L A I (4 7T r%) ] (1 - v r I c) , (2.48a-P-L)" 

Comment: 

If one would like to consider "rapidly rotating" particle, the result would be more complicated. So, 
we will consider that the particle does not rotate. The correct result for the flux of the incoming solar 
radiation interacting with the grain equals to 
Q' ext [ L Q A I (4 7T r%) ] (1 - 2 v r / c) , (2.48a-P-L-corrcct) . 

Q' ext = 2 for the special case treated by the authors. The second term in expression (2.48a-P-L-correct) 
accounts for the Doppler shift and the change of concentration of photons between the Sun's rest 
frame and that of the particle. 

The difference between Eqs. (2.48a-P-L) and (2.48a-P-L-correct) is evident. Incorrect form of Eq. 
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2. 

"The absorbed flux is reradiated isotropically and can be written as a mass loss rate in the particle's 

frame of motion (using E = mc 2 ): 

[ L A I (4 7T c 2 r%) ] (1 - v r I c) . (2.48b-P-L)" 

Comment: 

In general, the absorbed flux can be reradiated in an arbitrary way, not only isotropically. The mass 
loss rate, corresponding to the outgoing radiation in the particle's frame of motion, is: 
Q' ext [ L Q A I (4 7r c 2 r%) ] (1 - 2 v r / c) . (2.48b-P-L-correct) 

3. 

"As the particle moves relative to the Sun with velocity v, there is a momentum flux from the 
particle as seen in the rest frame of the Sun, since the particle emits more momentum in the forward 
direction than in the backward direction (Fig. 2.13-P-L). The momentum flux is equal to 
[ - L Q A I (4 7T c 2 r| ) ] (1 - v r I c) v , (2.48C-P-L) 

which can be generalized to the case in which the particle reflects and/or scatters some of the 
radiation impingent upon it via multiplication by Q pr ." 

Comment: 

As the particle moves relative to the Sun with velocity v, there is a momentum flux from the particle 

as seen in the rest frame of the Sun, since the particle emits more momentum in the forward direction 

than in the backward direction (Fig. 2.13-P-L). The momentum flux is equal to, for the case treated 

by Poynting (1903) and Robertson (1937), 

[ - 2 L e A I (4 7T c r|) ] [ (1 - v r / c) r - (1 / 2) x 

((1 - v r I c) r - (1 - 2 v r I c) v I c) ] , (2.48c-P-L-correct) 

see Eq. (76) in Klacka (2008b). The last result can be generalized to the case in which the particle 
reflects and/or scatters some of the radiation impingent upon it. However, this cannot be done in the 
simple way as stated by the authors (via multiplication by Q pr )- The correct result is given by Eq. 
(76) in Klacka (2008b). 



4. 
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"The net force on the particle in this more general case is given by: 

Frad = [ Lq Q pr A I (4 7T C t%) ] (1 - V r / c) f 

- [ L Q Q pr A v I (4 7T c 2 r%) ] (1 - v r / c) v (2.49a-P-L) 
^[L Q Q pr A/ (4vr c r 2 Q ) } [ (1 - 2 v r / c) r - (v e / c) § } . (2.49b-P-L) 

The first term in equation (2.49b-P-L) is that due to radiation pressure and the second and third 
terms (those involving the velocity of the particle) represent the Poynting-Robertson drag." 

Comment: 

The net force on the particle in this more general case is given by: 

F r ad = [ Lq Q ext A I (4 7T C T%) ] (1 - V r / c) T 

- [ L Q Q ext A I (4 7T c r|) ] [ (1 - v r / c) r - (Q pr / Q ext ) x 
((1 - v r I c) f - (1 - 2 v r I c) v I c) ] (2.49a-P-L-correct) 

~ [ L & Q pr A I (4tt cr%)}[(l -2v r / c)r - (v e / c) §] . 
The last force is the force due to the radiation pressure. 

6.27. Carroll and Ostlie (2007) 

We find the following definition in Carroll and Ostlie (2007, p. 806): 

"The Poynting-Robertson effect (a consequence of the headlight effect discussed in example 4.3.3) can 
cause ring particles to spiral in toward the planet. When particles in the rings absorb sunlight, they 
must re-radiate that energy again if they are to remain in thermal equilibrium. The original light was 
emitted from the Sun isotropically, but in the Sun's rest frame the re-radiated light is concentrated 
in the direction of motion of the particle. Since the re-radiated light carries away momentum as well 
as energy, the particle slows down and its orbit decays." 

Physics: 

The statement introduced above is one of the wide-spread ideas about the essence of P-R effect. 
Let us consider spherical particle. Re-emitted thermal radiation carries off energy and momentum per 
unit time (in the proper frame of particle) 



rpf si/ n f 

— ^abs^ ) 

P't = - F' e = , 



(54) 
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where C' abs is an absorption cross-section of the particle. I.e., the four- momentum of thermal radiation 
per unit time is 

p£ = \ C' abs S' (1 ; 0) . (55) 



Due to the emission of the energy, the rest mass of the particle will decrease: 

(56) 



fdm\ E' T C' abs S' 



T 



In the stationary frame, the four-momentum of thermal radiation per unit time will be 

P»t = \ C' abs S' ^ , (57) 

where is the particle's four-velocity. Thus, thermal radiation acts on the perfectly absorbing particle 
by four-force 

^r) = " Pt = ~ \ S' - . (58) 



dr ) T c aos c 

Using the fact that 

— I = I — u» + m [— — (59) 



dr J T \dr J T \ dr J T 

and also Eqs. (56) and (58), the four-acceleration of the particle is 

due to the thermal emission. Thus, the re-emitted thermal radiation produces no acceleration of the 
spherical particle. The particle is accelerated by incoming and by outgoing scattered radiation. But 
although the acceleration has the component in opposite direction to the particle's velocity v, the 
particle does not slow down, but its orbital velocity increases (see Sec. 6.21). 

6.28. Meyer et al. (2007) 

Meyer et al. (2007, p. 580) write: 

"The radial component of the radiation force is known as radiation pressure." ... "The tangential 
component of the radiation force is known as the Poynting-Robertson (P-R) drag. This acts on all 
grains and causes their orbits to decay into the star (where the grains evaporate) at a rate a = — 2 
a I a, where ..." 



Physics: 



In general, the radiation pressure force has three components, for arbitrarily shaped body. One 
radial and two nonradial. Equation of motion does not correspond to the P-R effect, in this general 
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case. The difference in motion of spherical and nonspherical dust grains was confirmed experimentally 
(Krauss and Wurm 2004). Only spherical grains behave in accordance with the P-R effect. 

The Poynting-Robertson effect is the effect of radiation pressure upon the motion of spherical 
particle. Nonradial components of the radiation pressure force equal zero, in this case. As for the 
P-R drag definition, much confusion exists in the literature, see Sec. 5. In any case, relativistically 
covariant equation of motion yields the P-R effect and its decomposition into several parts (classified 
as "effects") is nonphysical. 

6.29. Sykes (2007) 

Sykes (2007, p. 685) writes: "Comets were long thought to be the origin of the zodiacal cloud. 
However, estimates of dust production by short-period comets fell far short of that needed to main- 
tain the cloud in steady state against losses from particles spiralling into the Sun. This mechanism, 
where the absorption and reemission of solar radiation continually decrease particle velocity, is called 
Poynting-Robertson drag." 

Physics: 

What does physics say, in reality? Particle speed of the spiralling toward the Sun is an increasing 
function of time, according to Eqs. (46)-(48). The same states Eq. (45), since particle's semi-major axis 
a is a decreasing function of time and (3 is considered to be a constant. The same result follows from 
Eqs. (49)-(50): Eq. (50) yields decrease of eccentricity e and Eq. (49) produces an increase of the speed 
of motion for a decreasing eccentricity. Thus, the statement of the Sykes that "the absorption and 
reemission of solar radiation continually decreases particle velocity" does not correspond to reality. 
Moreover, totally reflecting particle is also described by the P-R effect - no absorption/reemission 
exists in this case. 

6.30. Km gel (2008) 

Kriigel (2008, p. 173-174) offers very similar explanation as Woolfson (2000), see also Sec. 6.9. He 
considers "a grain of mass to, radius a and geometrical cross-section a geo = na 2 circling a star at 
distance r with frequency u and velocity v = uor. The angular momentum of the grain is I — mr 2 uj 
= mrv" . Further he states:" Let M* and L* be the mass and luminosity of the star. Because there is 
mostly optical radiation, the absorption coefficient of the grain is C abs ~ a geo and the particle absorbs 
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per unit time the energy 
AE = L+a geo /(47rr 2 ) . 

In thermal balance, the same amount is reemitted. Seen from a non-rotating rest frame, the stellar 
photons that are absorbed travel in radial direction and carry no angular momentum, whereas the 
emitted photons do because they partake in the circular motion of the grain around the star. If we 
associate with the absorbed energy AE a mass m p hot — AE/c 2 , the angular momentum of the grain 
decreases per unit time through emission by 

dl/dt = - m phot rv = - (AE/c 2 ) rv = - L*<r geo l/ (4irc 2 r 2 m) . " 

Next procedure is in principle the same as the Woolfson's one. Kriigcl calculates a decrease of the 
particle's distance from the star based on the decrease of the orbital angular momentum of the grain. 

Kriigel also tries to explain the P-R effect from the view of the observer in a frame corotating with 
the particle (i.e. particle's proper frame). He states: "In such a reference frame, the stellar photons 
approach the grain not exactly along the radius vector from the star, but hit it slightly head-on in 
view of the aberration of light.... The photons thus decrease the angular momentum of the grain and 
force it to spiral into the star." 

Physics: 

From the view of the Sun. 
Calculation in Sec. 6.27 shows that the re-emitted thermal radiation from spherical grain acts on 
it by nonzero force. Thus, the re-emitted radiation changes the grain's orbital angular momentum. 
But this action is caused only by the decrease of the grain's mass due to the re-emission of 
the radiation and the corresponding acceleration equals zero. Thus, this phenomenon is not an 
essence of the Poynting-Robertson effect in the sense that electromagnetic radiation gives some 
acceleration to the particle. Moreover, thermal radiation does not represent the whole outgoing radia- 
tion which decreases particle's angular momentum. See Sees. 6.11 and 6.25 for more detail explanation. 

From the view of the particle. 
Kriigel's statement that the force originating from the incident radiation decreases (orbital) angular 
momentum of the particle is, similarly to Harwit's one (see Sec. 6.25), senseless. The velocity of 
the particle equals zero in its proper reference frame and thus the particle has no orbital angular 
momentum in this frame. 
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6.31. Liou and Kaufmann (2008) 

1. 

Liou and Kaufmann (2008, p. 426) write about "forces due to solar radiation and solar wind, including 
solar radiation pressure, PR drag, and solar wind drag." 

Comment: 

There are two physical forces: solar radiation force and solar wind force. These forces cannot be 
physically decomposed into various parts, as it is presented by the authors. This holds also for 
explanations presented in Liou and Kaufmann (2008, p. 426-427). Similarly, physically correct solar 
wind force should automatically contain possibility of nonradial velocity of the solar wind and time 
dependent properties of the solar wind, but no separation of these "effects" is physically acceptable. 

2. 

Liou and Kaufmann (2008, p. 426) write about the drag forces (PR drag and solar wind drag): 
"Although these drag forces are weak compared to the velocity-independent radial pressure force 
component, they dissipate energy and momentum and thereby cause the particles to eventually spiral 
into the Sun." 

Comment: 

The velocity-dependent terms of the radiation pressure and solar wind pressure forces dissipate total 
energy. But the momentum is not dissipated, since the speed of the particle increases during the 
process of spiralling toward the Sun. The magnitude of the momentum increases. The magnitude 
of the angular momentum decreases. Moreover, nonradial velocity component of the solar wind can 
increase particle's total energy and angular momentum (Klacka et al. 2008). 

3. 

Liou and Kaufmann (2008, pp. 426-427) write: "The PR drag can be thought of as arising from an 
aberration of the sunlight as ssen from the particle and a Doppler-shift induced change in momentum." 

Comment: 

As it is discussed in Sec. 7.4.5 in Klacka (2008b), the velocity-dependent term — v/c (v - heliocentric 
velocity of the particle, c is the speed of light in vacuum) in the P-R effect is not produced by the 
aberration of light. The same holds for the effect of solar wind. Moreover, the force for the solar wind 
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effect does not contain the term — v/u where u is the speed of solar wind particles (see Eq. 29 in 
Klacka and Saniga 1993). The effect of aberration is discussed also in Sec. 6.2. 

In reality, not only the Doppler effect causes the P-R effect. Also change of concentration of 
photons plays an equally important role (see also Klacka 1992a: Eqs. 31, 32, 36, 37, or, Sees. 2.4 and 
2.5; see also Klacka 1993b: Eqs. 2-5). 

4. 

Liou and Kaufmann (2008, p. 427) write: "Thus, forces due to collisions with solar wind particles are 
analogous to forces due to radiation." 

Comment: 

Yes, we can say that the force for the solar wind and solar electromagnetic radiation effects are 
analogous. However, one has to bear in mind that there are important differences. The force for 
the solar wind effect on dust particle contains decrease of mass of the dust particle (Klacka and 
Saniga 1993, Klacka 1993a, Kocifaj and Klacka 2008 - this physical difference is often not considered, 
although the change of mass of the particle is admitted; see, e.g., Dohnanyi 1978 - Sees. 2.1.3 and 
2.2.6, Mukai and Yamamoto 1982, Mann 2009 - Eqs. 7-10, 7-11), nonradial velocity of the solar wind 
(Klacka 1994, Klacka et al. 2008), time dependence of the solar wind (Svalgaard 1977). 

5. 

Liou and Kaufmann (2008, p. 427) write: "The contribution of the solar wind to the drag forces on 
dust particles is typically parametrized as a constant fraction sw of the PR drag." 

Comment: 

The statement of the authors is used in their equations: Eqs. (5-LK), (7-LK), (10-LK), (11-LK), (17- 
LK), (18-LK), (21-LK), (23-LK), i.e., the equations contain the term (3(1 + sw). The same access 
can be found, e.g., in Liou and Zook (1997), Holmes et al. (2003). However, this is not correct. The 
correct form is (except for the forces) (3(1 + sw/Q' pr ) as it is presented in Klacka (1994, Eq. 24 
although the value of sw is different from those used by Liou and Kaufmann), Klacka (2004: Eqs. 176, 
177, 182, 184, 186, ...), Klacka et al. (2008: Eqs. 4, 5, 15, 16). Moreover, the correct form for the force 
is even more complicated: 

Fdrag = - [ S AQ pr / (r§c) ] [ (1 + sw / Q pr ) (v-tq/c) f + (1 + x'sw/Q pr ) v/c] , 
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where the value of x' (> 1) depends on material properties of dust particle. The last equation is the 
correct equation, not Eq. (5) in Liou and Kaufmann (2008, p. 427). 



6.32. http://en.wikipedia.org/wiki/Poynting-Robertson_effect (last modification August 23, 
2008) 

In this very popular webpage we can find the following definition of Poynting- Robertson effect: 
"The Poynting-Robertson effect, also known as Poynting-Robertson drag, is a proces by which 
solar radiation causes a dust grain in the solar system to slowly spiral inward. The drag is essentially 
a component of radiation pressure tangential to the grain's motion." 



There is also an explanation of the effect: 
"From the perspective of the grain of dust circling the Sun the Sun's radiation appears to 
be coming from a slightly forward direction (aberration of light). Therefore the absorption of 
this radiation leads to a force with a component against the direction of movement. ... From the 
perspective of the solar system as a whole (...), the dust grain absorbs sunlight entirely in a radial 
direction, thus the grain's angular momentum remains unchanged. However, in absorbing photons, 
the dust acquires added mass via mass-energy equivalence. In order to conserve angular momentum 
(which is proportional to mass), the dust grain must drop into a lower orbit. 

Note that the re-emission of photons, which is isotropic in the frame of the grain, docs not affect 
the dust particle's orbital motion. However, in the frame of the solar system, the emission is beamed 
anisotropically, and hence the photons carry away angular momentum from the dust grain. It is 
somewhat counter-intuitive that angular momentum is lost while the orbital motion of the grain is 
unchanged, but this is an immediate consequence of the dust grain shedding mass during emission 
and that angular momentum is proportional to mass. 

The Poynting-Robertson drag can be understood as an effective force opposite the direction of 
the dust grain's orbital motion, leading to a drop in the grain's angular momentum. It should be 
mentioned that while the dust grain thus spirals slowly into the Sun, its orbital speed increases 
continuously. 

The Poynting-Robertson force is equal to: 
Fpr = Wv/c 2 = [r 2 /(4c 2 )]^GM s L2/ jR 5 

where W is the power of the incoming radiation, v is the grain's velocity, c is the speed of light, r 
the object's radius, G is the universal gravitational constant, M s the Sun's mass, L s is the solar 
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Comments: 

At first, we must point out that the Poynting-Robertson drag is not the same as the P-R effect, 
according to the standard approach used in literature (see also Sec. 5). The radial radiative pressure is 
also included into this effect. Therefore also the introduced equation for the P-R force is not correct, 
it corresponds only to the component proportional to —v/c, moreover not precisely (see Eqs. 21). 

We have to stress that the P-R effect is the effect describing motion of a spherical body (with 
spherically symmetric mass distribution) under the action of radiation pressure. It has no sense to 
give a special names to individual parts of the total force. At first, the confusion exists in scientific 
literature, and, moreover, the readers are utterly confused by such nonphysical steps. 

It has been already mentioned, for several times, that the aberration of the incident light does 
not lead to a force with component against the direction of particle's movement. For an explanation 
see Sec. 7.4.5 in Klacka (2008b) and also Sec. 6.2. 

From the author's explanation it follows that dynamics of the particle is influenced (its orbit is 
changed) only by the incident (absorbed; he does not consider complete interaction including also 
scattering) radiation because of conservation of particle's orbital angular momentum at absorption of 
this radiation. The outgoing (thermal re-emitted) radiation does not influence the particle's dynamics. 

The incident radiation really does not change particle's angular momentum (see Eqs. 8 and 18 and 
also Sec. 6.11) because the action of the increase of particle's mass is compensated by the action of 
deceleration of the particle. Thus, the author is partially correct in this sense. But his argumentation 
implies that the total acceleration of the particle should be proportional to Q' abs (particle's mass 
increase is proportional to this one) and not to Q' pr , as it actually is (see Eq. 22). Moreover, the 
statement of the author is not consistent with the fundamental condition represented by Eq. (16). 
Decisive role plays incoming and outgoing radiation. In reality, the incoming radiation generates 
particle's acceleration proportional to Q' ext = 2, for the case treated by the author (see Eq. 20). 

Equation of motion is the relevant physical result. The equation of motion contains (linear) mo- 
mentum vector of the particle, not its angular momentum vector. If we have equation of motion in 
disposal, it is not problem to obtain exact equations for the angular momentum of the particle. Really, 
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we can consider angular momentum of the particle 



r x Pi 



= m Hi 



where 




in 



is angular momentum per unit mass of the particle, and, similarly, 

L out = r x p ou t = m H ut , 

where 

Hout = r x v out ■ 

The results for these quantities can be summarized as follows: 
dL m /dt = , 

dH m /dt = - (\dm/dt\/m) H in , 



where \dm/dt\ — w 2 SA'Q' ext /c 2 (of course, v in — v ou t = v). These results immediately yield, e.g., 
that the magnitude of the vector Hi n decreases (particle spirals toward the central star) and the 
magnitude of the vector H ou t increases for the cases Q' pr IQ' ex t < 1 (the particle would spiral outward 
from the central star, if only "out" process would be important) - everything is consistent with the 
results obtained from linear momenta vectors (compare with Sec. 7.4.1 in Klacka 2008b). 

6.33. Zimmermann and Gurtler (2008) 

Zimmcrmann and Giirtlcr (2008, p. 307) write: "P-R-Effekt, das Phanomen, womach kleine sich urn 
die Sonne bewegende Teilchen infolge der Absorption von Sonnenstrahlung abgebremst werden, was 
zur allmahlichcn Annaherung an die Sonne fuhrt. Die von der Sonne emittierten Lichtquanten trager 
einen zu ihre. Energie proportionalen Impuls, der auf ein aborbierendes Teilchen iibertragen wird. 
Infolge der Bahngeschwindigkeit der umlaufenden Partikeln und der endlichen Lichtgeschwindigkeit 
ergibt sich bei der Absorption cin Aberrationseffckt. Der auf ein Teilchen iibertragene Impuls hat 
auBer einer von der Sonne wegge richteten radialcn Komponente eine der Bewegung des Teilchens 
entgegenge richtete tangentiale, abbrem sendwirkende Komponente. Jo klciner und massearmer eine 
Partikel ist, umso groBer ist des Verhaltnis des tangentialen Impulses gegeniiber dem Umlaufimpuls 
des Tielchens, umso schneller erfolgt die Abbremsung. Andererseits wachst mit atnehmendem 
Tielchenradius das Verhaltnis des radialen Strahlungsdrucks zur gravitativen Anzichurg durch die 
Sonne. Bei Tielchen groBer als etwa 0.1 /im iiberwiegt die Abbreusung, die Teilchen spieralen 
in die Sonne. Bei kleineren Partiklen iiberwiegt der Strahlungsdruck, sie werden von der Sonne 



dLout/dt = - (\dm/dt\/m) (Q' pr /Q' ext ) L ou t , 
dH ut/dt = + (\dm/dt\/m) (1 - Q' pr IQ' ext ) H, 



out 1 



weggetrieben." 
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Physics: 

f. 

The P-R effect holds for spherically symmetric bodies (see also laboratory measurements obtained 
by Krauss and Wurm, 2004). Solar electromagnetic radiation and gravity cause secular spiralling of 
spherical grains toward the Sun (if their radii are larger than w 0.1 fi m). 

2. 

The essence of the P-R effect is not an absorption of the incoming light. Even motion of the particle 
characterized with mirror-like spherical surface (specular reflection) is characterized by the P-R effect 
(see Sec. 7.3.3 in Klacka 2008b). 

3. 

The aberration of the incoming light does not lead to a force with component against the direction 
of particle's movement. For an explanation see Sec. 6.2 (or Sec. 7.4.5 in Klacka 2008b). 

4. 

The P-R effect is the effect of electromagnetic radiation pressure on moving spherically symmetric 
bodies. 

6.34. Mann (2009) 
1. 

Mann discusses "radiation pressure force" (Mann 2009, Sec. 7.3.1). Moreover, she discusses the 
"Poynting-Robertson effect", independently (Mann 2009, Sec. 7.3.3). Finally, she writes: "Particles in 
bound orbit about the star for which radiation pressure force is smaller than the stellar gravity force 
migrate toward the star due to the azimuthal component of the radiation pressure force and this is 
called the Poynting Robertson effect." (Mann 2009, Sec. 7.3.3). 

Mann (2009) makes difference between the radiation pressure force and the P-R effect. As we 
have already mentioned in Sec. 5, this access does not correspond to physics. The P-R effect is 
the radiation pressure force acting on moving (nonrotating) body with spherically symmetric mass 
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distribution. Moreover, the P-R effect contains not only the transversal velocity term, but also other 
terms, including radial velocity term. In reality, nonradial components of radiation pressure force 
are zero for spherical bodies. The nonradial components of radiation pressure force may exist for 
nonspherical particles (Klacka 2004, 2008b) - proper frame of reference of the particle is the relevant 
system for defining radial and nonradial components. 

2. 

Mann (2009, p. 201) states: 

"The stellar wind force in the frame of the moving particles depends on the dust velocity v as: 
F sw =F sw [(l-v r/{rv sw )) r/r - v/v sw ] , (7.10-M) 

where F sw is the force on the dust for v — and v sw is the bulk velocity of the wind. The non-radial 
term in (7.10) is referred to as plasma or pseudo Poynting- Robertson drag force." (We have corrected 
the evident errorneous term in the original equation, Eq. (7-10-M) is correct.) 

Physics: 

i) 

The solar- wind force acting on a body is not of the form presented by Mann (2009 - see Eq. 7.10-M). 
The correct form is given by Klacka and Saniga (1993 - Eq. 29): 
Fsw =F sw [(l-v r/(rv sw )) r/r - x' v/v sw } , (29-K-S) 

where x' denotes the part of the incident energy which is lost from the grain during the interaction 
with the solar wind (measured in the reference frame of the grain; 1 < x' < 3, approximately, 
depending on material properties of the particle). The force is even more complicated if nonradial 
component of the solar wind velocity is included. 

ii) 

What is the sense of the words "plasma or pseudo Poynting-Robertson drag force"? We have just 
shown that the last term of the real force Fsw does not correspond to that presented by Mann 
(2009 - Eq. 7.10). Thus, the term "plasma or..." is misleading. And what about the rest velocity term? 

3. 

Mann (2009, p. 202) states: 

"The ratio of plasma P-R drags force to (photon) P-R drag force can be written as 



Klacka et al.: Poynting- Robertson effect 59 

(F sw /F rad ) c/v sw « ... (7.1 1-M) 

The factor c/v sw results from the difference of the aberration angles for photons or solar wind 
particles." 

Physics: 

We are aware that this kind of explanation is presented also in Burns et al. (1979, p. 12), or, Leinert 
and Griin (1990, pp. 226-227) - see also Sec. 6.5. However, the solar-wind force does not produce the 
term — F sw v/v sw (see Eq. 29-K-S presented above). Thus, explanation of the term as a consequence 
of the aberration is incorrect. Even the light-term — F p h v/c is not produced by the aberration of 
light, see Sec. 6.2 and also Sec. 7.4.5 in Klacka (2008b). 

6.35. Minato et al. (2004) 

We have already mentioned an analogy between the solar electromagnetic and corpuscular radiation 
in Sees. 6.30, 6.34. We will present detail comments on the paper by Minato et al. (2004), now: the 
incorrect statements on the solar wind effect may enable better understanding of the P-R effect. 

1. 

Minato et al. (2004 - Sec. 2.1) state that the contribution of solar-wind forces to the motion of dust 
grain is described as 

F sw [ (1 - 2 r/v sw ) e R - (r9/v sw ) e e } , (1-M) 

where v is the velocity of the grain v = r e# + r9 eg and v sw is the heliocentric solar-wind speed. 
Physics: 

The solar-wind force acting on a body is not of the form presented by Minato et al. (2004 - Sec. 
2.1, sec Eq. 1-M). The correct form is given by Klacka and Saniga (1993 - Eq. 29): 
F sw [ (1 - (1 + x 1 ) r/v sw ) e R - (x' r9/v sw ) e g } , (29-K-S) 

where x' denotes the part of the incident energy which is lost from the grain during the interaction 
with the solar wind (measured in the reference frame of the grain; x 1 > 1). The force is even more 
complicated if one uses also the nonradial component of the solar wind velocity. 



2. 
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"The second term on the right-hand side of Eq.(l-M) expresses the pseudo PR drag which causes 
dust grains to spiral into the Sun." (Minato et al. 2004 - Sec. 3.2) 



Physics: 

i) 

There does not exist the term — r8/v sw eg in the solar wind force. The real term is multiplied by the 
value x'. This is one of the differences between the P-R effect and the solar wind effect. Thus, there 
does not exist a simple analogy to the " P-R drag" . 
ii) 

There is another inconsistency in defining the "P-R drag". We have already discussed (see Sec. 5) 
that various definitions are used. Now, Minato et al. (2004 - Sec. 3.2) present another definition for 
the "P-R drag". Thus, according to the authors, the "P-R drag" is neither the term — F p h [ 2 (r/c) 
e R + (r0/c) eg ], nor the term — F p h [ (r/c) en + (rO/c) eg ], but the term — F p h (rO/c) eg. And 
what about the rest velocity terms in the P-R effect? How are they called in the " physics" of Wyatt 
(2009), Mann (2009) and Minato et al. (2004)? 

Really, the radiation pressure force for a body with spherically distributed mass corresponds to the 
P-R effect. Relativistically covariant formulation of the equation of motion shows that any separation 
of the equation into various (ad- hoc) components is of no physical sense and the term "P-R drag" 
should be omitted in future. 

By the way, the term — F sw (r0/v sw ) eg may not cause that dust grain spirals into the Sun. Also 
the outspiralling from the Sun may occur if the solar wind bombardment significantly decreases mass 
of the grain. Moreover, non-radial solar wind velocity may be more important than the term — F sw 
(rd/v sw ) eg. 

3. 

"The ratio of 7 of the pseudo P-R drag to the P-R drag is given by 7 = (F sw /F ph ) c/v sw , where the 
factor c/v sw comes from the difference of the aberration angles for solar wind and light." (Minato et 
al. 2004 - Sec. 3.2) 



Physics: 



Klacka et al.: Poynting- Robertson effect 



61 



The solar- wind force does not produce the term — F sw r9/v sw eg. Thus, explanation of it as a 
consequence of the aberration is incorrect. Even the light-term — F ph (r6/c) e g is not produced by 
the aberration of light, see Sec. 6.2 (Sec. 7.4.5 in Klacka 2008b). 

7. Simple explanation of the P-R effect 

Astronomers often want to present a simple explanation of a physical effect. We have just presented 
many incorrect explanations of the P-R effect. The correct explanation can be found in Klacka (2008b). 
However, very simple explanation can be very useful for scientists and students. 

7.1. Simple explanation without equations 

Physical correct and simple explanation with no equation can be as follows: 

Poynting- Robertson effect - The effect of electromagnetic radiation pressure on a moving spherical 
particle. The incident radiation can be represented by a parallel beam of photons and the particle is 
a nonrotating body of spherically symmetric mass distribution. 

Let the Sun (or another star) be the only source of radiative and gravitational forces acting on 
the particle. The particle of radius smaller than about 0.1 micron (the value depends on the particle's 
optical properties and the luminosity of the star) is expelled from the Solar System (the star) while 
the larger particle slowly spirals toward the Sun (the star). 

Solar radiation pressure force acting on small dust particles (less than 0.1 micron in radii) is larger 
than the solar gravitational force. As for larger particles, the conservation of particle's mass, during 
particle's interaction with the incoming solar radiation, is important. The conservation of the mass 
produces the term proportional to — v/c in the radiation pressure force acting on the particle; v is 
particle's heliocentric orbital velocity and c is the speed of light. The drag/frictional term proportional 
to — v/c causes decrease of particle's total energy (and angular momentum) and the particle slowly 
spirals toward the Sun. The magnitude of the momentum (and the speed of the particle) slowly 
increases. 

7.2. Mathematical description 

Homogeneous beam of parallel photons interacts with non-rotating spherical body/particle. If the 
incoming momentum per unit time is p' 'incoming in the proper frame of reference of the particle, then 
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the total outgoing momentum per unit time is 

P outgoing (l ^pr/^ext) P incoming • 

C' pr and C' ext are pressure and extinction cross-sections. The cross-sections are calculated from Mie's so- 
lution of Maxwell's equations (Mie 1908). The following inequalities hold: C' pr > and < C' pr /C' ext < 
2. The case of "perfectly absorbing spherical particle" (treated by Poynting 1903 and Robertson 1937) 
is, in reality, characterized by C pr — nR 2 , C' ext — 2irR 2 , where R is radius of the particle. The force 
acting on the particle is 

dp / (It p incoming P outgoing — {^pr I ext) P incoming ; 

where r is proper time measured in the reference frame of the particle. Using the fact that 

P incoming — i,S C ex ^/c) €. , 

where S 1 is the flux density of radiation energy (energy flow through unit area perpendicular to the 
ray per unit time), c is the speed of light and e' is unit vector describing direction and orientation of 
the travelling photons with respect to the particle. The last two equations yield 
dp'/dr = (S'C' pr /c) e' . 

This equation represents an electromagnetic (light) pressure acting on the particle. Moreover, the mass 
of the particle does not change and this corresponds to the condition for the energy of the particle E' 
dE'/dr = . 

The left-hand sides of the last two equations are components of the energy-momentum four-vector 
(four-momentum). This enables us to find the equation of motion of the particle in any frame of 
reference, if we use Lorentz transformation. If the particle moves with velocity v in a frame of reference, 
then the equation of motion of the particle is 
dp^/dr = (w 2 SC pr /c) (&" - u?/c) , (*) 

where = (E/c ; p) is four-momentum of the particle, w = 7 (1 — v ■ e/c), 6 M = (1/w ; e/w), 
= (7c ; "fv) is four-velocity of the particle, e is unit vector describing direction and orientation of the 
travelling photons with respect to the stationary frame of reference (laboratory frame) in which the 
particle moves with velocity v, and, finally, 7 = l/-\/l — (v/c) 2 is the Lorentz factor. 

The equation of motion (*) is known as the Poynting-Robertson effect. The Poynting-Robertson 
effect is the effect corresponding to the electromagnetic (light) pressure. 



To the first order in v/c, Eq. (*) reduces to 
dv/dt = [SC' pr /(mc)] [ (1 - v ■ e) e - v/c } . 



Remark: 
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The process of interaction between the incident radiation and the spherical particle can be treated 
as a two-step process. 



The incoming radiation acts on the particle with a force 

(d P »/dT) in = (w 2 SC' ext /c) , 

and the force causes the following acceleration of the particle 
{du»/dr) m = [ w 2 SC' ext /(mc) } (&" - u"/c) , 
or, to the first order in v/c: 

(dv/dt) in = [ SC' ext /(mc) ] x [ (1 - v ■ e/c) e - v/c] . 



The outgoing radiation acts on the particle with a force 
(dp»/dT) out = - (w 2 SC' ext /c) [ (C' pr /C' ext ) u"/c + (1 - C' w /C' ext ) b» } , 
and the force causes the following acceleration of the particle 
{du»/dT) out = - [ w 2 SC' ext /{mc) ] (1 - C" pr /C" ext ) (6" - u"/c) , 
or, to the first order in v/c: 

{dv/dt) out = - [ SC' ext /{mc) } (1 - C' pr /C' ext ) x[(l-v e/c) e - v/c] . 



It can be easily seen that (dv/dt) out oc + v if C'pr/C'ext < 1. This holds also for the case treated 
by Poynting and Robertson {C' pr /C' ext — 1/2). 

Often the dimensionless efficiency factors of extinction and pressure are used: Q' ext — C' ext /{nR 2 ), 
Q' pr — C' pr /(irR 2 ), where R is radius of the particle. 

Orbital evolution of interplanetary spherical dust particle under the action of solar electromagnetic 
radiation and the solar gravity is thoroughly understood and explained in Klacka (2004) and Klacka et 
al. (2007). Orbital evolution of spherical particle under the action of solar gravity and electromagnetic 
radiation leads to a secular decrease of semi-major axis and eccentricity. Shift of perihelion may exist 
already in the case when the equation of motion is considered to the first order in v/c, if optical 
properties of the particle change. 
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8. Secular evolution of orbital elements 

Let us consider orbital evolution of a body under the action of gravity and electromagnetic radiation 
of a central star (e. g., Sun). The action of electromagnetic radiation is taken to be the P-R effect. 
Accuracy to the first order in v/c is considered. 

We have already discussed that physics does not admit to divide the P-R effect into partial com- 
ponents. Thus, we have to take into account that the P-R effect represents the disturbing force to the 
Keplerian motion given by the gravity of the central star. We are interested in secular evolution of 
orbital elements (short-term oscillations in orbital elements are not important for us). We will proceed 
in accordance with Klacka (1994) and Klacka (2004 - Sec. 6). 

If parameter (3 changes, then only numerical calculations can be realized. If equation of motion in 

the form 

GM 
v = — r 

„ GM ( / v ■ e\ v i 

is numerically solved, then orbital elements can be calculated on the basis of Eqs. (47) in Klacka (2004 
- the correct right-hand side of the last equation contains v ■ ey/(e \J GM Q Jp ) — 1/e, if M = M Q , 
instead of ... — 1). 

In any case, an approximation of constant (3 can be useful, at least for comparison. This case can be 

treated in a fully analytical approach, if secular evolution of orbital elements is searched. The results 

can be summarized as follows, as for semi-major axis a and eccentricity e: 

-l 

— dx , 



1 f**h + 0(l + 4 + 2e () ca8x)/(l-e%) 

a = ap(l-4)' — / ± ^ — ^ 

^ Jo 1 + eg COS £ 



3/2^ [ 2 * ^(l-P) 2 e} + ^-2p(l-(3)e cosx 
' ' i? 2tt J (l + e/3 cosa;) 2 



(1 + e/3 cosxY 

) eg cosx 

dx , (62) 



where 

dap_ _ _ pGM 2 + 3e| 



dt ' C 9 \ 3 / 2 

ag 



4, =-2 13 — 77, (63) 



(i-4) 

deg 5 a GM eg 

al (l - e^ 

Orbital elements with subscript (3 correspond to the central Keplerian acceleration — GM(l — (3)r/r 3 
- non- velocity radial term of acceleration caused by radiation force in Eq. (61) is added to the gravi- 
tational acceleration of the star (however, it is often stated that Eqs. 63 hold for central acceleration 
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- GMr/r 3 , see, e. g., Burns et al. 1979 - Eqs. 30a, 30b, 30c, 30d, 31a, 31b, 46a, 46b, 47, 48, Mignard 
1992). The set of differential equations (Eqs. 63) is the set presented by Robertson (1937), Wyatt and 
Whipple (1950). Moreover, initial conditions for Eqs. (63) are required. If the particle is released from 
the parent body (quantities relating to the parent body are denoted by subscript P) with velocity A 

A = Av R e PR + Av T e PT + Av N e PN , (64) 

while 

v = vr epR + vt epT , 



vr = JGMpp 1 ep sin(6p — up) 



v T = \jcM pp 1 [1 + e P cos(e P -lop)} ; (65) 
Vp = a P (l — e p)> Pp = a /3 (l — e p^j an d epN = epii x epr, then (compare Gajdosfk and Klacka 



1999; Klacka 2004 - Eq. 59): 



Pp 



pp (v T + Av T ) + (Av N ) 



2 



1-/3 v 2 



T 



2 x 2 



2 / 1 + e P cos/p v TS 

Gf-i in — I ^ ^ 9 t 



\ 1-/3 v 2 



1 + e P cos/p v TS \ f v R + Avr Y 
1-/3 v T J V v t J 

Pj3 in 



1 - e 2 
1 e /3 in 



«T + Avt . Avn „ . . 

COS Iff i n = COSlp — cosBp smi P , 

VTS VTS 

v T + Av T . . n 

Sim 8 i n COS i 1/3 in = smip COSiZp + 

VTS 

^ VN (cosQp cos ip cos fip — sin9p sin Op) 

VTS 

v T + Av T . . . n 

Sim fi in Smilg in = sin Ip SlXlilp + 

VTS 

~" ' (cos0p cos ip sin Op + sin Op cos ftp) 



VTS 

ep in COS } 'pin = (1 + e P cos/p) - 1 



ep m sin fp 



Pp 

1 + ep cos/p vts (vr + &vr) 
1-/3 4 



wt + Av T . . „ 
sinifl j„ cosfc>/3 j„ = simp cos Hp + 

COSlp , 

sinz^ ^ sin Qp in = sinip sinOp , 
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v% s = (v T + Av T ) 2 + (Av N ) 2 . 



(66) 



The set of equations represented by Eqs. (62)-(66) fully corresponds to detailed numerical calcu- 
lations of vectorial equation of motion, if we are interested in secular evolution of eccentricity and 
semi-major axis for the case when central acceleration is defined by gravity alone. We have to stress 
that "ep i n < 1" and "ep does not correspond to pseudo-circular orbit" are tacitly assumed. 




Fig. 1. Dependence of a/ ap as a function of secular value of ep. Dashed lines depict functions (1 — ep) 
I [1 — ep + (3(1 + ep)] and (1 + ep) / [1 + ep + (3(1 — ep)\. The dashed lines confine the zone where 
instantaneous ratio of semi-major axis (central acceleration — GMr/r 3 ) and ap can occur. 

It is worth mentioning that instantaneous time derivatives of semi-major axis and eccentricity 
defined by central acceleration — GMr/r 3 may be both positive and negative, while secular evolution 
yields that semi-major axis a and eccentricity e are decreasing functions of time. As for the central 
acceleration — GM(l—(3)r/r 3 , the instantaneous time derivative of semi-major axis is always negative, 
while the instantaneous time derivative of eccentricity may be both negative (near apocentrum) and 
positive (secular evolution yields that semi-major axis ap and eccentricity ep are decreasing functions 
of time - see Eqs. 63). 
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Fig. 2. Dependence of e as a function of secular value of ep. Dashed lines depict functions |(1 — (3)ep — 
f3\ and (1 — (3)ep + (i. The dashed lines confine the zone where instantaneous eccentricity (central 
acceleration — GMr/r 3 ) can occur. 



For the special case A = Eqs. (66) reduce to 



1 + ep cos/p 



ap i n = a P (1 - /?) 1 - 2/3 



6/3 ir. 



/ 1 1-4 -2/3 (l + e P cos/p) 



(i-/?r 

where fp = 6p — u)p, up ,„ has to be obtained from 

(i + e P cos f P 

e/3 in COS (Bp - m ) = — , 

ep m sin (Bp - wp in ) = ^ - - , 



(67) 
(68) 



p in 



^P , ip in=ip , 0/3 



e 



p • 



(69) 
(70) 



Physics of Eqs. (67)-(70) is following: meteoroid escapes from the Solar System when the orbital energy 
becomes positive and this can happen when the energy due to the radial component of the radiation 
force is included, without it being necessary for this force to exceed the gravitational attraction (Harwit 
1963). Some figures may be found in Klacka (1992b, 1993c). As an example we may mention that 
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time 



Fig. 3. Time evolution of semi-major axes a and ap. T is the time of spiralling toward the star. The 
value of CLi n is given by the first of Eqs. (62). The following time limits hold: lim^T a = lim t ^T ap — 
0. 



particle of = (1 — e P )/2 moves in parabolic orbit if ejected at perihelion of the parent body, and, in 
an orbit with eccentricity ep in = |1 — 3ep|/(l + e P ) if released at aphelion; ep in = for = e P = 
1/3 and aphelion ejection. 

We may mention that a little more simple procedure is obtained when semi-major axis dp is re- 
placed by pp = ap(l — ep) (a quantity referred to as semi-latus rectum). Then pp = pp in {ep/ep i«) 4 ^ 5 , 
and: 

3/2 

(71) 



1 - e l in {PP/P0 in) 



dpp = _ 2p GM 
dt c pp 



de _ 5 GM e 8 J 5 in (l - ' 

dt 9 r n 2 3/5 

ai L c Pp in eJ 



5/2 



3/2 



(72) 



Eqs. (62)-(63) show that secular evolutions of semi-major axes a, ap and eccentricities e, ep are 
decreasing functions of time. This means that eccentric orbits become more circular, during the parti- 
cle's orbital evolution (constant is assumed). This corresponds to the fact that apocentric distances 
(aphelion distances) Q, Qp decrease more rapidly than the pericentric distances (perihelion distances) 
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Fig. 4. Time evolution of eccentricities e and ep. T is the time of spiralling toward the star. The 
following time limits hold: lim^T e = /3, linit_>T ep = 0. 

q, qp, as it is evident also from the equations 



dx 



1 - yf{\-ff)'e\ + (3 2 - 2(3 (1-/3) ep cos x 



1 — + (3 (l + e 2 + 2ep cos x^j (1 + ep cos x) 



(73) 



1 f 

Q = ap (l - e^) 5 2 — y Z Q , 

1 + ^(1 - (if el + f3 2 - 2(3 (1 - J) ep cosx 



1 — + (3 (l + e 2 + 2ep cos x^j (1 + e/3 cos x) 2 
where Eqs. (63)-(66) can be used (or, Eqs. 71-72, too), or 

a P = \ + Qp) , 

= 1 - qp/Qp 
1 + qp/Qp ' 
and qp, Qp are given by the following equations: 

= - (3 — J£_ 1 b 1P + 3 Q/3 
dt c 2Qp y/qpQp qp + Qp 



(74) 



(75) 
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dQp = _ pGM_ Qp_ 1 3gp + 5Qp 
dt c 2qp ^JqpQp qp + Qp 

dq^JdQA = (M_Y( 1 _ 2 Qe-v ) < l (76) 



dt ) ' \ dt ) \Qp) \ Zqp + bQp 

For completeness we have to introduce initial conditions 

Pp in 



1(5 in — 



1 + e/3 i n 
Pp in 



Qp in = Z , (77) 

J- - ep in 

where pp in and ep i n are given by Eqs. (66). 

It is worth mentioning that instantaneous time derivatives of pericentric and apocentric distances 
defined by central acceleration — GMr/r 3 may be both positive and negative, while secular evolution 
yields that q and Q are decreasing functions of time. As for the central acceleration — GM(1 — {3)r/r 3 , 
the instantaneous time derivatives of pericentric and apocentric distances are both less than or equal 
zero (secular evolution yields that qp and Qp are decreasing functions of time - see Eqs. 76). 

The time of spiralling of the particle into the star is 

6.4 x 10 2 M sun ("g m [AU]f(l-e} m ) p in z 3/5 

T[yrs] = -M L • (78) 

(see also Eq. 51). 

Figs. 1 and 2 depict the behaviour of Eqs. (62). The dashed lines correspond to the dependences 
(1 — ep) j [1 — ep + (3(1 + ep)] and (1 + ep) j [1 + ep + (3(1 — ep)] for the ratio of instantaneous values 
of semi-major axes and (1 — (3)ep — (3\ and (1 — (3)ep + (3 for instantaneous values of eccentricities 
corresponding to the central acceleration — GMr/r 3 (see Eqs. 96-97 in Klacka 2004). Time evolution 
of semi-major axes and eccentricities is presented in Figs. 3 and 4. The time of spiralling toward the 
central star is given by Eq. (78). While the values of a, ap and ep approach to zero, the value of e 
tends to the value (3 (see Eqs. 62-63). 



9. Conclusion 

We have presented physical comments to " explanations" of the Poynting- Robertson effect published 
during the past thirty years by almost hundred of scientists working in the field. We have used mainly 
scientific papers and monographs. Although we have only restricted access to literature, we hope 
that physics of other published presentations can be covered by our paper. Moreover, we hope that 
our paper, together with papers by Klacka (1992a, 1993b, 2004, 2008a, 2008b), can help in better 
understanding the dynamical interaction of electromagnetic radiation with a moving body. 
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Moreover, we have added two simple explanations of the P-R effect (see Sec. 7). They are not only 
simple, but they concentrate the physical essence of the phenomenon. These simple explanations (or 
their modifications) may be used in scientific literature, encyclopediae and textbooks. 

The confusion in the scientific literature about the term called the " Poynting- Robertson drag" 
(see, e.g., Dohnanyi 1978, Gustafson 1994, Minato et al. 2004, Wyatt 2009, Mann 2009) can be easily 
removed if one takes into account that the relativistically covariant form of the P-R effect does not 
enable any division into partial (covariant) terms. The correct physics yields the P-R effect as a whole 
effect of the radiation force on a moving body and no "P-R drag" exists as a partial relativistically 
covariant equation of motion. 

The difference between the effects of solar electromagnetic and corpuscular (solar wind) radiation 
is stressed (see Sees. 6.34, 6.35). 

Sec. 8 presents secular evolution of orbital elements of a spherical particle if the action of gravity 
of a central star and it's electromagnetic radiation is considered. 

The physical understanding of the P-R effect enables correct derivation of the action of stellar winds 
on evolution of dust grains. This can help in physical explanation of the observed dust distribution at 
stars with planetary systems. 
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